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y  This  thesis  refines  an  earlier  demonstrated 


methodology  which  combines  response  surface  methodology  and 
experimental  design  concepts  to  describe  the  output  of  a 
deterministic  model.  In  an  effort  to  maximize  the  combat 
power  delivered  to  a  theater  commander,  several  critical 
factors  are  varied  to  determine  their  impact  upon  the 
delivery  of  men  and  equipment  to  the  objective  area.  The 
effects  of  time,  materials  handling  equipment  ''t'Miffi-) — ^ 
capability,  airfield  ramp  space,  distance  between  the 
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Aerial  Port  of  Debarkation  'TAFOB-)^  and  Forward  Operating 
Location  fFOht^the  availability  and  capabilities  of 
different  aircraft,  and  the  impact  of  aircraft  attrition 
are  examined  in  the  model.  This  model  does  not  use  metric 
ton-miles  per  day  as  a  measure  of  effectiveness,  but  rather 
it  places  a  time  dependent  value  on  the  delivery  of  combat 
units  to  the  objective  area.  The  deterministic  model 
includes  168  equations  and  288  different  variables.  A 
response  surface  was  generated  using  an  appropriate 
experimental  design,  with  the  response  surface  predictions 
averaging  greater  than  a  96  percent  level  of  accuracy.  The 
analysis  section  demonstrates  the  advantages  and 

\ 

disadvantages  of  the  sensitivity  analysis  normally 
associated  with  linear  programming,  and  the  advantages  and 
limitations  of  using  response  surface  methodology  to 
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generate  a  multidimensional  analysis  of  the  model.  While 
the  optimum  airlift  force  mix  may  be  heavily  scenario 
dependent,  this  research  provides  additional  insight  into 
the  complexities  of  developing  an  appropriate  airlift  force 
to  meet  future  Department  of  Defense  requirements. 


THE  USE  OF  MATHEMATICAL  PROGRAMMING  AND  RESPONSE 


SURFACE  METHODOLOGY  IN  OPTIMIZING  THE  AIRLIFT  FORCE 
STRUCTURE  FOR  A  FAR  EASTERN  THEATER  OF  OPERATION 


I.  Introduction 


Background 

In  1975,  the  commander  of  the  Military  Airlift 
Command  stated  the  enhancement  of  our  airlift  capability  is  a 
critical  step  in  acquiring  an  efficient  airlift  force. 
(5:2-10)  In  order  to  protect  its  vital  security  interests, 
the  United  States  must  be  able  to  support  its  forces 
overseas.  In  the  early  stages  of  a  major  conflict,  airlift 
will  play  a  dominant  and  vital  role  in  delivering  the  men  and 
equipment  required  for  the  operation.  An  effective  airlift 
force  must  be  able  to  meet  strict  time  constraints  and  it 
should  be  able  to  deliver  combat  units  directly  from  their 
place  of  origin  to  the  employment  area.  The  current  airlift 
force  structure  does  not  possess  that  capability  (4:iii). 

This  type  of  concern  precipitated  the  Congressionally 

Mandated  Mobility  Study  (CMMS),  which  documented  the 
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requirement  for  additional  airlift  capability,  especially  of 
outsized  cargo.  As  a  short  term  solution,  in  January  of  1982 
the  Department  of  Defense  announced  the  purchase  of  50 
improved  C-5B  aircraft.  When  the  last  C-5B  is  delivered  in 
1989,  long  range  intertheater  airlift  capability  will  reach 
46  million  ton-miles  per  day,  still  26  per  cent  short  of  the 
66  million  ton-miles  per  day  goal  established  by  the  CMMS 
(24:28-29) . 

As  an  additional  factor,  the  remaining  airlift  fleet  is 
in  need  of  modernization.  The  C-141,  the  current  workhorse 
of  intertheater  airlift,  is  over  20  years  old  and  will  either 
need  to  be  replaced  around  the  turn  of  the  century  or  be 
upgraded  with  a  life  extension  program.  The  primary 
intratheater  aircraft,  the  C-130,  has  several  models  which 
are  rapidly  approaching  the  end  of  their  service  life 
(21:14).  Airlift  planners  must  now  determine  how  to 
effectively  meet  the  need  for  a  replacement  aircraft  capable 
of  fulfilling  the  vital  intratheater  mission. 

In  1980,  the  report  of  the  Secretary  of  the  Air  Force  and 
Air  Force  Chief  of  Staff  to  Congress  stated  that  "the  current 
airlift  deficiency  is  judged  to  be  the  greatest  problem  the 
Air  Force  faces  in  executing  the  national  military  strategy" . 
(34:17) 

While  the  CMMS  was  still  in  progress,  airlift  planners  at 
the  Air  Staff  and  the  Military  Airlift  Command  formed  a  C-X 
Task  Force  to  develop  the  operational  concept  and  to 
determine  what  capabilities  would  be  desirable  in  a  future 


aircraft.  This  Task  Force  produced  the  Preliminary  System 
Operational  Concept  (PSOC)  which  described  the  intended 
purpose  of  the  C-X.  The  PSOC  developed  several  scenarios  and 
focused  upon  the  mission  required  of  the  aircraft  for  each 
scenario.  The  mission  requirements  were  determined  and 
passed  to  the  aircraft  manufacturers ,  who  were  given  the 
flexibility  to  design  the  aircraft  to  meet  the  needs  of  the 
Air  Force.  Boeing,  Lockheed,  and  McDonnell  Douglas  all 
responded  to  the  request  for  proposal,  and  in  1981,  the 
Source  Selection  Evaluation  Board  convened  to  examine  all  the 
proposals.  In  August  of  the  same  year  the  Secretary  of  the 
Air  Force  selected  the  McDonnell  Douglas  design  as  the 
winner.  The  design  has  subsequently  been  designated  the 
C-17  (24:28). 

The  C-17’ s  ability  to  deliver  cargo  from  onload  bases  in 
the  United  States  directly  to  the  operational  forces  in  the 
field  adds  a  new  dimension  to  airlift  responsiveness  and 
highlights  why  acquisition  of  the  C-17  has  become  a  major 
priority  for  MAC  (1:37).  The  flexibility  of  the  C-17  in 
meeting  the  intertheater,  intratheater,  and  outsize  cargo 
carrying  requirement  has  made  the  C-17  the  cornerstone  of 
MAC’S  effort  to  meet  future  airlift  requirements.  MAC  must 
now  convince  Congress  that  the  C-17  is  a  better  choice  than 
the  C-5B,  even  though  50  C-5B’s  have  already  been  approved  by 
Congress . 

The  major  disadvantage  to  the  C-17  is  its  acquisition 

cost.  At  over  $100  million  per  aircraft,  the  added 
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capability  does  not  come  cheap.  The  C-17  can  deliver 
outsized  cargo  directly  to  the  objective  area,  but 
battlefield  commanders  might  be  reluctant  to  commit  such  an 
expensive  resource  to  a  hostile  environment  (1:51).  While 
the  C-5B  costs  over  $140  million  per  aircraft,  it  can  carry 
twice  as  much  cargo  by  volume,  36  pallets  versus  18  for  the 
C-17,  and  40  per  cent  more  cargo  by  weight,  242,000  pounds 


K 


B 


versus  172,000  pounds  for  the  C-17.  Yet,  in  a  report  from 
the  GAO  to  the  Chairman  of  the  Committee  on  Armed  Services, 
the  GAO  stated  that  the  Air  Force  feels  that  neither  the  C-5A 
or  the  C-5B  can  routinely  and  safely  land  or  takeoff  from 
small  austere  airfields.  Although  Lockheed  disagrees,  the 
Air  Force  feels  that  these  conditions  would  be  operating  the 
crew  and  aircraft  near  their  limits  with  small  margin  for 
error.  This  report  stressed  the  improved  capability  of  the 
C-17,  such  as  its  smaller  size,  better  maneuverability, 
better  combat  offload  capability,  and  airdrop  capabilities, 
and  how  these  improvements  are  vital  to  the  accomplishment  of 
the  entire  airlift  mission  (34:2).  Due  to  its  size  and 
design  characteristics  the  C-5B  would  be  restricted  to 
operating  within  only  well  developed  airfields  and  would 
probably  not  operate  in  a  hostile  environment.  Outsized 
cargo  would  have  to  travel  overland  to  get  to  the  objective 
area  which  means  less  threat  to  the  aircraft  but  a  longer 
time  delay  involved  in  getting  the  required  equipment  to  the 
battlefield  commanders.  As  stated  by  Colonel  Bullard,  during 


his  assignment  to  the  Air  War  College,  "unless  the  men, 


equipment,  and  supplies  are  delivered  to  the  combat  area, 
only  a  marginal  utility  will  be  achieved  by  transporting  them 
to  the  general  theater  of  operation."  (4:23)  There 
definitely  exists  a  tradeoff  between  capability,  costs, 
effectiveness,  and  vulnerability,  and  the  Air  Force  and 
Congress  must  agree  upon  which  factors  have  priority. 

Research  on  this  topic  and  methods  capable  of  quantifying 
the  relative  advantages  of  each  aircraft  have  been  discussed 
only  peripherally. 

After  reviewing  such  factors  as  military  utility,  program 
costs,  manpower  requirements,  force  stabilization,  and  force 
modernization,  the  CMMS  recommended  a  force  structure  relying 
heavily  upon  the  C-17.  The  CMMS  stated  this  airlift  force 
would  provide  the  United  States  with  the  best  ability  to 
rapidly  deploy  and  sustain  its  fighting  forces  in  pursuit  of 
national  objectives  (8:V-12). 

A  paper  submitted  at  the  U.S.  Army  War  College  assessed 
the  capabilities  of  the  C-17  had  it  been  available  for  use  in 
two  actual  airlift  operations.  The  first  operation,  AHUAS 
TARA83,  was  a  combined  operation  in  Honduras,  and  the  second, 
Urgent  Fury,  was  the  U.S.  rescue  mission  to  Granada.  In  both 
cases  the  C-17  was  shown  to  be  more  effective  and  efficient 
than  the  aircraft  actually  used  in  the  operations  (13:14-29). 

In  a  paper  submitted  at  the  Air  Command  and  Staff 
College,  several  advantages  and  disadvantages  of  the  C-17 
were  discussed.  In  spite  of  reservations  about  using  the 

C-17  in  a  moderate  threat  environment,  this  paper  recommends 
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purchasing  the  maximum  number  of  aircraft  (approximately  300) 
(1:56-57) . 

As  further  evidence  of  the  utility  of  the  C-17,  a  thesis 
submitted  by  Major  John  Stone  stressed  the  "direct  delivery" 
concept  of  the  C-17.  In  his  conclusion  he  stated  that  the 
opponents  who  denounce  the  C-17  because  its  payload  capacity 
is  less  than  the  C-5  do  not  understand  that  cargo  capacity  is 
of  little  value  if  the  cargo  is  not  delivered  in  a  timely  and 
efficient  manner  to  where  it  is  needed  (32:91). 

As  pointed  out  by  Major  Soligan,  while  assigned  to  the 
Air  Command  and  Staff  College,  airlifts’  role  in  the  Airland 
Battle  still  needs  to  be  researched  and  defined.  He 
emphasized  the  importance  of  the  Airlift  Concepts  and 
Requirements  Agency  (ACRA)  and  its  role  in  meeting  with  the 
theater  commanders  to  resolve  the  proper  role  of  airlift  in 
the  deployment  and  resupply  effort  required  to  develop 
doctrine  and  support  the  desired  concept  of  operations 
(30:53) . 

None  of  the  articles  thus  far  have  dealt  with  the  issue 
of  airlift  force  structure  using  an  analytical  or  simulation 
model.  In  1984  Army  Captain  James  Cook,  an  AFIT/GST  student, 
demonstrated  his  methodology  using  Sequential  Goal 
Programming  and  Response  Surface  Methodology  to  determine  how 
the  interaction  of  unit  weight,  combat  attributes,  logistics 
needs,  and  airlift  resources  could  be  jointly  optimized  for  a 
given  contingency.  By  developing  a  model  which  he 

demonstrated  using  212  variables  and  136  constraints,  he  was 
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able  to  analytically  determine  the  optimal  force  mix  or 
incremental  advantage  in  deployed  power  attainable  by  an 
incremental  change  in  airlift  resources  (7:4). 

It  is  the  purpose  of  this  investigation  to  expand  and 
refine  the  methodology  used  by  Captain  Cook  to  optimize  the 
airlift  force  structure  required  in  a  Far  Eastern  theater  of 
operation. 

Problem  Statement 

There  is  currently  no  analytical  model  which  accurately 
and  fully  quantifies  the  airlift  force  mix  which  would 
optimize  the  combat  power  delivered  to  a  battlefield 
commander  in  a  Far  Eastern  theater  of  operation. 

Research  Objective 

The  primary  objective  of  this  research  effort  is  to 
refine  an  earlier  demonstrated  methodology  that  quantifies 
the  value  of  C-5  and  C-17  aircraft  in  intertheater  airlift 
and  the  C-130  and  C-17  in  an  intratheater  role.  The 
capability  of  each  aircraft  to  meet  mission  requirements  is  a 
function  of  several  major  factors  including  the  quantity  and 
capabilities  of  the  specific  aircraft  and  the  survivability 
of  the  aircraft  in  a  threat  environment.  The  impact  of  the 
different  factor  levels,  including  quantity  of  aircraft 
available,  quantity  of  materials  handling  equipment  (MHE) 


available,  distance  between  APOD  and  FOL,  available  ramp 
space,  and  aircraft  attrition,  upon  the  measure  of 
effectiveness  of  the  model  and  the  relative  value  of  the  C-5, 
C-130,  and  C-17  needs  examination  and  quantif ication  in  order 
to  determine  an  efficient  and  effective  airlift  force  mix  for 
a  specific  scenario. 

S.ggJgg 

The  determination  of  an  appropriate  airlift  force  mix  is 
heavily  influenced  by  the  scenario  requirements.  The 
selection  of  a  Far  Eastern  theater  of  operation  is  intended 
to  represent  the  type  of  scenario  which  the  United  States 
must  be  prepared  to  meet;  a  scenario  which  combines  such 
critical  factors  as  intertheater  and  intratheater  distances, 
airfield  and  MHE  saturation,  different  aircraft  capabilities, 
and  the  presence  of  a  threat  environment.  Several 
limitations  to  the  scenario  are  listed  in  the  scenario 
development  chapter.  The  linear  programming  model  and  the 
response  surface  equations  are  intended  to  demonstrate  the 
methodology  of  effectively  utilizing  RSM  and  experimental 
design  concepts  to  explain  the  output  of  a  deterministic 
model.  While  necessarily  restrictive,  the  methodology  is 
intended  to  provide  valuable  insight  into  the  complexities  of 
determining  airlift  force  structure  requirements. 


Earns  t 


The  format  for  this  thesis  consists  of  the  following 
major  areas.  Chapter  II  covers  the  general  methodology  by 
reviewing  response  surface  methodology  and  experimental 
design,  followed  by  a  discussion  of  recent  applications  of 
these  concepts.  Chapter  III  develops  the  formulation  of  the 
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168  constraints  and  288  variables  required  for  the  linear 
programming  model.  The  next  chapter,  Chapter  IV,  introduces 
the  data  base  used  to  develop  the  scenario.  Chapter  V  is  an 
analysis  of  the  results,  including  the  development  of  the 
experimental  design  and  response  surface  equations.  Chapter 
VI  contains  a  summary  of  conclusions  and  a  list  of  possible 
recommendations  for  further  research. 


The  overall  objective  of  this  study  is  to  refine  the 


methodology  established  by  Captain  Cooke  in  order  to 
maximize  the  combat  power  delivered  to  a  Far  Eastern 
theater  of  operation  as  a  function  of  time.  While 
considering  such  factors  as  aircraft  capabilities, 
scenario  limitations,  and  attrition  of  aircraft,  this 
objective  can  be  stated  as  one  overall  goal:  maximize 
combat  power  delivered  while  minimizing  the  time  required 
to  deliver  the  army’s  assets  to  the  objective  area. 
Response  Surface  Methodology 

Response  Surface  Methodology  (RSM)  is  a  mathematical 
method  which  permits  a  particular  set  of  mathematical  and 
statistical  tools  to  be  used  by  researchers  to  solve 
certain  types  of  problems,  especially  in  industrial 
research.  This  concept  has  been  used  extensively  in 
industry  to  examine  factor  relationships  and  to  determine 
which  combinations  of  factors  result  in  the  optimum  output 
(29:58). 

Mathematical  models  formulated  by  researchers  consist 
of  a  number  of  input  variables,  or  independent  variables, 
which  describe  the  output  or  response  of  a  given 
experiment.  The  response  surface  generated  by  the 
mathematical  equations  provides  a  picture  of  the  influence 


of  the  various  input  variables  upon  the  primary  measure  of 
interest  (response  variable). 

Ideally,  the  researcher’s  model  would  include  all  the 
different  variables  that  affect  the  output  response. 
Realistically,  constructing  an  appropriate  mathematical 
model  to  represent  a  real  world  problem  will  most  probably 
be  an  iterative  process  where  variations  to  the  model  will 
eventually  result  in  the  model  response  approximating  the 
value  of  the  real  response.  In  a  complicated  process,  the 
potential  for  not  including  a  critical  factor  in  the 
development  of  a  deterministic  model  can  be  significant. 
The  exclusion  of  one  or  more  critical  factors  in  a  model 
produces  a  model  with  biased  results.  The  factor  or 
factors  not  considered  in  the  model  bias  the  model  results 
and  produce  conclusions  that  differ  from  actual  results. 

A  good  mathematical  model  will  tend  to  minimize  bias  error 
as  well  as  the  second  type  of  modeling  error,  variance 
error. 

Variance  error  is  caused  by  an  inability  to  hold 
certain  factors  constant  during  an  experiment.  In  the 
chemical  industry,  the  need  to  minimize  variance  error  is 
critical  to  the  successful  interpretation  of  a  given 
experiment.  The  method  of  least  squares  regression  allows 
the  researcher  to  minimize  error  by  determining  whether 
error  is  due  to  an  incomplete  list  of  the  factors  or  to 
poor  measurement  methods . 


RSM  is  used  here  to  determine  the  least  number  of 
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measurements  which  are  needed  to  accurately  fit  the  model 
to  the  data.  If  RSM  can  be  used  to  accurately  determine 
an  appropriate  mathematical  model  to  represent  the  output 
of  an  experiment,  then  the  researcher  is  no  longer 
required  to  run  the  entire  experiment  in  order  to  analyze 
possible  outcomes  from  new  combinations  of  factors.  As 
long  as  factor  combinations  are  chosen  which  remain  within 
the  bounds  of  the  experiment,  the  response  of  the 
experiment  can  be  closely  approximated  using  the 
mathematical  model. 

Box  and  Wilson  first  developed  the  theory  of  Response 
Surface  Methodology  in  1951.  Throughout  the  years  several 
useful  experimental  designs  have  been  developed  to  expand 
the  usefulness  of  RSM  for  various  purposes.  The  most 
common  and  easily  understood  designs  are  for  first-order 
and  second-order  models.  A  first-order  model  is  linear  in 
form.  It  is  said  to  be  linear  because  no  parameter 
appears  as  an  exponent  or  is  multiplied  or  divided  by 
another  parameter.  The  independent  variable  is  also 
considered  linear  because  it  appears  only  in  the  first 
power.  Second-order  models  are  slightly  more  complex  and 
include  the  use  of  quadratic  variables  and  variables 
representing  two-way  interactions.  The  following 
equations  illustrate  first  and  second-order  models: 

Y  =  value  of  response  variable 

Bq  and  B^  =  parameters  (variables) 
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x  =  value  of  independent  variable 
i  =  random  error  term 
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In  analyzing  the  coefficients  of  a  given  model,  it  is 
desirable  to  have  coefficients  which  are  uncorrelated. 
Correlation  refers  to  the  influence  one  variable  has  upon 
the  value  of  another  variable  in  the  equation.  In  1960 
Box  and  Behnken  suggested  possible  new  designs  which 
possessed  a  high  degree  of  orthogonality .  Orthogonal 
variables  are  uncorrelated;  therefore,  the  various  terms 
in  the  mathematical  equation  are  independent  of  each 
other.  This  allows  the  various  coefficients  to  be 
compared  to  each  other  independently  or  as  ratios.  In  the 
Box  and  Behnken  designs  only  the  constant  term  Bq  and  the 
quadratic  estimates  are  correlated  to  one  another,  but 
current  data  has  shown  that  the  correlation  between  these 
coefficients  has  been  very  small  and  has  not  presented  any 
major  problems  in  practical  applications  (29:60). 

As  stated  previously,  the  error  from  a  model  is 
produced  by  a  combination  of  bias  error  and  variance 
error.  Variance  error  reflects  a  measurement  error 


whereas  bias  error  is  caused  by  an  improperly  fitting 
model  probably  caused  by  the  omission  of  critical 


variables  in  the  mathematical  equations  describing  the 
experiment . 

In  the  past,  most  of  the  discussion  on  experimental 
design  has  centered  around  the  importance  of  minimizing 
this  variance  error.  Very  little  work  has  been 
accomplished  in  the  area  of  minimizing  bias  error.  Myers 
summarizes  the  work  accomplished  by  Box  and  Draper  in  his 
book  "Response  Surface  Methodology”.  (22) 

The  earlier  applications  of  RSM,  which  emphasized 
variance  reduction,  centered  around  chemical  experiments 
and  the  need  to  minimize  the  error  caused  by  improper 
methods  of  measurements  in  the  laboratory.  The  chemical 
industry  has  made  good  use  of  RSM  in  predicting  outcomes 
and  in  determining  the  effects  of  interaction  upon  output 
responses . 

As  the  use  of  RSM  increases,  there  is  a  growing  need 
to  emphasize  the  importance  of  minimizing  bias  error. 
Analytical  math  models  which  use  the  same  inputs  produce 
results  with  no  real  variance.  The  outputs  of  these 
deterministic  models  do  not  change  if  the  input  values  do 
not  change;  therefore,  variance  is  not  a  factor.  Bias 
error  now  becomes  the  prevalent  type  of  modeling  error, 
and  research  continues  into  methods  of  minimizing  it. 
Another  important  factor  to  consider  when  constructing  a 
response  surface  is  selecting  an  appropriate  experimental 
design. 


When  selecting  an  experimental  design  to  construct  a 


response  surface,  several  guidelines  need  to  be 
emphasized.  As  stated  by  Myers,  the  factors  must  be 
quantitative  and  continuous,  the  function  must  be  able  to 
be  approximated  using  a  low-order  polynomial,  and  the 
independent  variables  x^,  • • •*£  should  be  highly 

controlled  and  be  measurable  with  only  a  negligible  amount 
of  error  (22:62). 

Other  requirements  for  a  good  experimental  design 
include  the  ability  to  accurately  estimate  the 
coefficients  while  not  requiring  an  inordinate  number  of 
observations.  In  order  to  properly  fit  a  second-order 
model,  the  variables  in  the  model  must  be  evaluated  at 
three  different  levels,  usually  at  a  low,  high,  and 
intermediate  value  (20:462).  Another  important  feature  is 
the  designs’  rotatability .  As  stated  by  Myers,  "A  design 
is  said  to  be  rotatable  when  the  variance  of  the  estimated 
response,  that  is,  the  variance  of  "y" ,  is  a  function  only 
of  the  distance  from  the  center  of  the  design  and  not  on 
the  direction.”  (22:139)  As  a  result,  points  in  the 
factor  space  which  are  the  same  distance  from  the  origin, 
are  treated  as  being  equally  important.  Or  stated  another 
way,  rotatability  insures  a  symmetric  generation  of 
information  in  the  space  defined  by  the  variables 
(3:457). 


Another  desirable  feature  of  an  experimental  design 


is  to  produce  coefficients  which  are  uncorrelated.  As 
stated  in  the  section  on  response  surfaces,  uncorrelated 
coefficients  are  produced  by  orthogonal  designs.  In  an 
orthogonal  design  the  columns  of  the  X  matrix  are 
orthogonal  to  each  other  so  in  a  completely  orthogonal 
design  the  (X’X)  matrix  would  produce  zeroes  for  all  the 
off-diagonal  elements.  The  coefficients  would  therefore 
be  independent  and  uncorrelated  with  each  other.  The 
three-level  rotatable  central  composite  design  selected 
for  this  thesis  was  developed  by  Box  and  Behnken  (3:461) 
and  is  nearly  orthogonal  as  evidenced  by  the  correlation 


matrix  listed  in  Table  X.  Only  the  constant  term  and 


the  quadratic  terms  are  correlated  with  each  other  and  as 
previously  stated,  the  degree  of  correlation  is  extremely 
small  and  does  not  affect  the  significance  of  the  design 
output  (29:60).  The  coded  and  non-coded  designs  are 
listed  as  Appendix  A  and  Appendix  B  respectively.  The 
next  section  of  this  chapter  discusses  several  recent 
applications  which  have  incorporated  the  concepts  of 
response  surface  methodology  and  experimental  design. 
Military  Applications  of  RSM  and  Experimental  Design 


Within  the  last  several  years,  a  number  of  papers 
have  been  written  which  applied  the  methods  of  RSM  and 
Experimental  Design  to  various  military  problems.  The 
next  few  pages  will  discuss  the  implications  and 
conclusions  of  several  of  these  efforts. 
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One  of  the  first  important  applications  which 
combined  RSM,  mathematical  programming,  and  experimental 
design,  was  written  for  SAC  in  1975  by  LTC  Smith  and  Dr. 
Joseph  Mellichamp.  They  established  a  methodology  which 
will  be  employed  in  this  thesis,  which  allows  for 
multidimensional  impact  analysis  of  military  problems. 

The  article,  which  was  first  published  in  1975,  pointed 
out  the  classical  sensitivity  analysis  of  varying  one 
parameter  at  a  time  has  some  major  disadvantages.  First, 
only  one  parameter  can  be  evaluated  at  a  time.  Second, 
the  current  method  does  not  allow  for  any  numerical 
ranking  of  the  importance  of  each  parameter  to  the  total 
solution.  Third,  no  information  can  be  provided  about  the 
interrelationships  that  may  exist  between  the  different 
factors  under  study .( 29 : 56 )  RSM  will  allow  a 
multidimensional  analysis  of  the  same  problem.  Using  the 
correct  experimental  design  to  optimally  select  a  number 
of  combinations  of  the  input  parameter  values,  linear  or 
quadratic  programming  can  find  the  optimum  value  of  the 
response  variables.  Linear  regression  can  then  be  used  to 
determine  the  coefficients  of  the  mathematical  equation 
defining  the  response  surface  (29:58).  As  long  as  the 
factor  combinations  are  within  the  range  of  the  parameters 
used  to  develop  the  initial  equation,  the  mathematical 
equation  can  then  be  used  as  a  method  of  providing 
sensitivity  analysis  for  the  true  response  surface. 


The  second  paper  to  be  discussed  is  the  thesis  of 
Captain  James  Cooke.  His  work  forms  the  foundation  for 
this  thesis  effort.  He  developed  the  methodology  which 


will  be  refined  and  expanded  to  determine  the  optimum 
combat  power  delivered  to  a  specific  theater  of  operation, 
for  a  specified  time  period,  and  within  certain 
sustainability  levels.  He  related  the  effects  of 
intertheater  and  intratheater  movement  and  the  deployable 
unit's  capabilities  to  the  optimization  of  combat  power 
delivered  to  a  given  objective  area.  Unit  weights,  combat 
power  attributes,  resources,  and  in-theater  constraints 
were  used  to  find  global  optimum  solutions  (7:4). 

Response  Surface  Methodology  was  used  as  a  method  of 
sensitivity  analysis  by  examining  the  effects  of  several 
independent  variables  upon  combat  power  delivered  to  the 
front.  The  flexible  response  surface  was  generated  using 
a  fractional  factorial  design  based  upon  the  need  to 
maintain  orthogonality  between  variables  and  to  minimize 
bias  error  within  the  model.  Captain  Cooke  demonstrated 
his  methodology  by  using  an  RSM  model  containing  five 
variables,  which  consisted  of  the  number  of  C-5A  and  C-17 
aircraft,  APOD  size,  materials  handling  equipment,  and 
distance  from  the  APOD  to  the  forward  operating  location. 
He  analyzed  the  incremental  advantage  in  deployed  combat 
power  attainable  by  an  incremental  change  in  airlift 
resources  (7:4).  It  is  the  purpose  of  this  thesis  to 
apply  Cooke’s  methodology  to  a  specific  scenario  In  an 
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effort  to  further  examine  the  impact  of  certain  variables 
within  the  mathematical  model,  the  additional  variables  of 
C-130  availability  and  aircraft  attrition,  and  the  impact 
of  aircraft  attrition  upon  the  amount  of  combat  power 
delivered,  will  be  examined  in  the  response  surface 
model.  Previous  discussions  in  airlift  force  structuring 
have  indicated  the  importance  of  the  scenario  in 
determining  the  optimal  force  mix,  and  this  scenario  will 
center  around  the  deployment  of  combat  units  from  CONUS  to 
the  Republic  of  Korea. 

The  different  potential  applications  of  RSM  can  be 
demonstrated  by  another  thesis  presented  in  December  of 
1984  by  Captain  Manacapilli.  He  combined  the  use  of  RSM, 
economic  production  functions,  economic  theory, 
deterministic  models,  and  lagrangian  techniques  to  help 
identify  cost  effective  strategic  force  mixes.  His 
methodology  was  used  to  fit  economic  production  functions 
to  the  response  surface  of  an  unclassified  linear 
programming  nuclear  exchange  model.  The  thesis  presented 
a  methodology  which  combined  the  applications  of  economics 
with  RSM  to  analyse  the  output  of  a  deterministic  model. 
His  research  showed  that  if  economic  production  functions 
are  properly  fitted  to  a  response  surface,  the  information 
gained  can  provide  useful  insight  into  the  production 
process  itself  (15:6.4). 


Another  example  of  the  far  reaching  impact  of 
combining  mathematical  programming  with  response  surface 
methodology  was  the  work  done  by  Captain  Graney,  who 
examined  the  use  of  RSM  as  an  aid  to  multi -dimensional 
decision  making.  Graney’ s  study  presents  the  optimization 
methodology  for  a  deterministic,  mathematically  modeled 
system  as  defined  by  the  consecutive  levels  of  one  measure 
of  effectiveness  (MOE)  and  the  optimization  of  a  second 
goal  or  mission.  This  methodology  allows  the  decision 
maker  the  option  of  a  large  number  of  alternatives .  It  is 
capable  of  comparing  the  MOE’s  throughout  the  entire  range 
of  alternatives  and  allows  the  use  of  non-commensurate 
measures  of  effectiveness  (11:79).  Graney  demonstrates 
his  methodology  by  examining  an  aggregated  military  force 
structure  problem  as  defined  by  a  linear  programming 
arsenal  exchange  model.  The  two  measures  of  effectiveness 
are  counterforce  and  countervalue,  and  the  system  consists 
of  five  weapons  systems  and  ten  target  classes.  As  noted 
by  the  author,  the  number  of  feasible  alternatives  for  the 
system  exceeds  three  billion,  but  optimal  effectiveness  of 
these  alternatives  can  be  approximated  through  the  use  of 
response  surface  methodology. 

Using  an  appropriate  fractional  factorial  design,  a 
response  surface  was  generated  for  each  MOE.  The 
methodology  allows  for  a  tradeoff  between  the  two  MOE’s  by 
optimizing  the  second  measure  of  effectiveness  as  it  is 
constrained  to  the  contour  lines  of  the  first  measure  of 


effectiveness.  This  thesis  demonstrated  the  usefulness  of 
RSM  in  examining  the  alternatives  in  a  multi-criteria 
comparison  of  alternatives  and  approximating  the  optimal 
effectiveness  of  those  alternatives. 

Three  recent  thesis  efforts  have  attempted  to  develop 
interactive  user-friendly  computer  models  to  assist  users 
in  applying  and  analyzing  complex  deterministic  or 
simulation  models. 

The  first  model  was  developed  by  Captain  Tate  in 
response  to  Captain  Cooke's  thesis  recommendation  that  his 
mathematical  model  be  developed  into  a  user-friendly 
computer  package. 

Captain  Tate  developed  an  interactive  goal 
programming  model  designed  to  analyze  the  rapid  deployment 
of  units  of  the  armed  forces  to  a  selected  crisis  area. 

The  mathematical  model  used  by  Tate  in  creating  his 
computer  model,  DEPLOY,  was  taken  largely  from  the 
mathematical  model  developed  by  Captain  Cooke,  discussed 
earlier  in  the  text.  Tate  developed  some  alternative 
methods  for  describing  Cooke’s  mathematical  equations,  but 
for  the  most  part,  his  contributions  center  around  his 
creation  of  a  user-friendly  model  which  would  capture  the 
complex  mathematical  relationships  established  by  Cooke. 
Tne  allocation  routine,  PAGP,  needs  refinement  in  order  to 
correctly  interface  with  the  computer  package  DEPLOY,  and 
this  enhancement  was  presented  in  the  thesis  as  a 
recommendation  for  further  research  (33:7-4). 


Nevertheless,  Tate  did  develop  Cooke’s  mathematical  model 
into  a  computer  package  which  analyzed  the  accomplishment 
of  Army  objectives  as  well  as  determined  the  minimum  force 
structure  required  to  accomplish  the  goals  set  by  the 
Army.  A  potential  advantage  of  this  model  is  its  ability 
to  perform  two  types  of  analysis.  The  first  determines 
the  degree  to  which  the  Army’s  goals  can  be  met,  given 
that  the  number  and  types  of  available  aircraft  and 
deployable  units  are  predetermined  parameters .  The  second 
type  of  analysis  determines  the  minimum  force  structure 
required  to  meet  the  combat  power  goals  set  by  the  Army 
(33:7-2) . 

The  second  computer  model,  developed  by  Mr.  Meitzler, 
implemented  a  user-friendly  program  module  into  a 
relatively  new  deterministic  simulation  model, 
Budget/Readiness  Analysis  Technique  (BRAT).  The  model  was 
developed  by  the  Studies  and  Analysis  Directorate  of  the 
Acquisition  Logistics  Center  (AFALC/XRS),  and  it  converts 
the  logistics  resources  of  maintenance  manpower,  support 
equipment,  and  spares  into  a  sortie  rate  to  determine  a 
system's  state  of  mission  readiness.  He  applied  the 
techniques  of  response  surface  methodology,  experimental 
design,  and  multiple  regression  to  develop  a  module  which 
uses  low  order  polynomials  to  approximate  the  response 
surface  of  the  BRAT  model .  The  response  surface  can  be 
used  to  provide  the  capability  to  quickly  perform 
sensitivity  analysis  on  the  logistics  resources  and 
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simulation  parameters  (18:vii,24).  He  discusses  the 
importance  of  selecting  the  correct  experimental  design 
based  upon  any  prior  knowledge  of  the  form  of  the  response 
surface.  First  and  second  order  designs  were  examined, 
and  the  importance  of  the  design  defining  the  boundaries 
of  the  experimental  region  was  examined.  A  first  order 
surface  would  use  a  2  design  where  each  k-factor  is  taken 
at  two  levels,  and  a  quadratic  or  second  order  design 
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using  a  3  design  would  insure  each  K-factor  is  taken  at 
three  levels.  The  major  contribution  of  this  thesis  lies 
in  the  implementation  of  a  user-friendly  computer  module 
into  the  BRAT  module  which  provides  the  user  the 
capability  to  use  RSM  to  analyze  the  outputs  of  the  model 
and  to  perform  sensitivity  analysis  on  the  logistics 
resources  and  simulation  parameters  (18:24). 

A  third  interactive  computer  package  was  developed  by 
lLt  Sparrow.  She  designed  a  user-friendly  response 
surface  methodology  computer  package  which  could  be 
attached  to  any  Fortran  based  simulation  model.  The  model 
was  designed  to  produce  a  response  function  to  accurately 
predict  the  relationship  between  the  input  variables  and 
the  response  variable.  The  computer  package  was  developed 
to  be  implemented  on  the  VAX  11/780,  which  can  run  the 
simulation  models  and  any  required  Fortran  subroutines. 

The  package  was  designed  to  be  user-friendly  and  offered 
five  different  design  types,  but  it  did  not  include  the 
Box  and  Behnken  designs  which  will  be  used  for  this  thesis 
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(31:106).  Their  inclusion  into  the  program  was  listed  as 
a  possible  enhancement  to  the  current  model . 

The  last  several  papers  discussed  the  capabilities  of 
response  surface  methodology  to  predict  the  output  of  a 
complex  model.  In  order  for  RSM  to  achieve  this  goal,  the 
response  surface  must  be  developed  with  the  appropriate 
experimental  design.  The  importance  of  selecting  an 
appropriate  design,  containing  certain  important 
characteristics,  was  analyzed  by  Captain  Ishihara. 

He  evaluated  the  effectiveness  and  efficiency  of 
several  bias  and  variance  minimizing  experimental 
designs.  He  used  several  experimental  designs  to  develop 
response  equations  for  a  deterministic  nuclear  exchange 
model.  The  following  criteria  was  used  to  evaluate  the 
design  points:  1)  the  total  number  of  design  points,  2) 
the  number  of  terms  in  the  response  equation,  3)  the 
accuracy  of  the  fit  of  the  response  equation,  4)  the 
orthogonality  of  the  design,  and  5)  the  rotatability  of 
the  design  (14: ix).  He  also  evaluated  the  various 
equations’  ability  to  estimate  the  true  surface  and  to 
predict  future  responses  from  a  given  set  of  input  data. 
This  paper  concluded  that  a  second  order  polynomial  can  be 
an  excellent  choice  for  estimating  the  response  surfaces 
of  deterministic  nuclear  models.  The  level  of  accuracy  of 
a  modified  equation  depends  upon  the  number  of  design 
points  and  upon  the  degree  of  rotatability  of  the  design 
(14:4-2).  The  response  equation  should  contain  as  few 


cross  product  and  squared  terms  as  possible.  The 
experimental  design’s  accuracy  should  improve  as  the 
number  of  design  points  increased,  as  long  as  the  design 
points  were  chosen  from  throughout  the  entire  surface. 
Additionally,  if  the  region  being  investigated  is 
spherical,  a  rotatable  design  improved  the  accuracy  of  the 
overall  response  equation  (14:4-2).  Another  conclusion  of 
Ishihara  is  that  specific  conditions  must  be  met  in  order 
for  the  coefficients  of  a  modified  response  equation  to 
accurately  represent  the  marginal  effects  of  change  of  the 
dependent  variable  caused  by  incremental  changes  to  the 
independent  variables.  These  conditions  include:  1)  a 
high  degree  of  multicol linearity  should  not  exist  between 
the  independent  variables,  2)  the  main  effect  terms  of  the 
response  equation  should  dominate  the  squared  an$i 
interaction  terms,  and  3)  the  modified  response  surface 
should  be  able  to  accurately  represent  the  true  surface 
(14:4-3).  These  conclusions  reaffirm  the  importance  of 
incorporating  a  good  experimental  design  into  the 
development  of  the  response  surface. 

All  of  these  papers  demonstrate  the  potential 
analytical  value  of  combining  response  surface  methodology 
with  an  appropriate  experimental  design  to  develop  an 
accurate  and  efficient  method  of  examining  the  effects  of 
critical  factors  upon  the  outputs  of  deterministic  models. 


This  chapter  will  be  used  to  develop  and  explain  the 
mathematical  expressions  describing  the  goals  and 
constraints  used  in  the  model.  The  formulations  of  the 
goals  and  constraints  for  this  model  are  taken  largely 
from  the  thesis  of  Captain  Cooke.  The  variables  and 
notation  for  considering  the  constraints  over  a  period  of 
days,  rather  than  for  every  day,  are  taken  from  the 
previous  work  of  Captain  Tate.  While  Cooke  formatted  his 
constraints  for  a  daily  operation,  Tate  transformed  those 
constraints  to  be  applicable  over  a  given  period  length. 
There  are  minor  differences  for  a  variety  of  constraints 
and  the  data  base  describing  military  hardware 
capabilities,  and  these  differences  will  be  fully  examined 
in  this  chapter.  The  variables,  multipliers,  and 
subscripts  will  be  discussed  first,  followed  by  an 
examination  of  the  development  of  the  goals  and 
constraints.  After  each  constraint  is  derived,  the  number 
of  that  type  constraint  required  for  the  model  will  be 
given,  followed  by  the  formulation  as  it  applies  to  this 
particular  model.  The  end  of  the  chapter  will  include  a 
summary  of  all  the  goals  and  constraints  required  for  this 
particular  model. 


The  variables  as  formulated  by  Captains  Cooke  (7:46- 
76)  and  Tate  ( 33 : chap3 , 1-31 )  are  defined  as  follows: 
Variables  Used: 

P(k,l):  variables  used  to  signify  excess  units 

at  the  FOL  which  will  be  available  for  use  at 
the  start  of  the  next  period.  "k"  refers  to  the 
cargo  type  and  “1"  refers  to  the  period. 

PS(k,l):  variables  used  to  signify  excess  supplies 
at  the  FOL  which  will  be  available  for  use  at 
the  start  of  the  next  period. 

A(k,l):  variables  used  to  signify  excess  units  at 

the  APOD  available  for  use  at  the  start  of  the 
next  period. 

AS(k,l):  variables  used  to  signify  excess  supplies 
at  the  APOD  which  will  be  available  for  use  at 
the  beginning  of  the  next  period. 

SI:  the  amount  of  supplies  transshipped  in  period  1 
from  the  APOD  to  the  FOL.  [ tons/period] 

x(i,j,k,l):  the  number  of  delivered  aircraft  loads, 
[ac /period] 

u(y,m,l):  the  number  of  units  deployed . [units/period] 

ATT(i):  the  expected  overall  attrition  for  type  "i" 
aircraft 
Subscripts  Used: 

a:  APOD 

i:  a  specific  aircraft  type,  range  is  from  i  to  I . 


For  this  model  i  ranges  from  1  to  6  with  the 
following  aircraft: 

1=  C-5B  Galaxy 
2=  C-17 

3=  C-141B  Starlifter 

4=  CRAF  Boeing  747 (cargo  type) 

5=  CRAF  McDonnell  Douglas  DC-8 (personnel  type) 
6=  C-130D/E  Hercules 

j:  a  specific  type  of  airlift  mission,  range  is  from 
1  to  J.  For  this  model  j  ranges  from  1  to  4. 

1=  Delivery  to  the  APOD  (Intertheater) 

2=  Direct  delivery  to  the  front  (Intertheater) 
3=  Airdrop  missions  (Intertheater) 

4=  Intratheater  missions 

k:  a  specific  cargo  type,  range  is  from  1  to  K.  For 
this  model  k  ranges  from  1  to  4. 

1=  outsize 
2=  oversize 

3=  bulk( including  supplies) 

4=  personnel 

1:  a  specific  period,  range  is  from  1  to  L.  For  this 
model  1  ranges  from  1  to  4 ,  with  each  period 
covering  5  days  for  20  days  total, 
m:  the  mode  of  delivery,  range  is  from  1  to  M.  For 
this  model  m  ranges  from  1  to  3. 

1=  delivery  to  the  front 

2=  delivery  to  the  APOD  and  remain  at  the  APOD 
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or  move  to  the  front  on  their  own 
3  =  delivery  to  the  APOD  and  move  to  the  front 
via  intratheater  airlift. 

y:  a  specific  unit  type,  range  is  from  1  to  Y.  For 
this  model  range  is  from  1  to  8 . 

1=  Airborne  units  of  the  82nd  Airborne  Division 

2=  HQ  units  from  the  82nd 

3=  Air  Assault  Units (AH-64  equipped) 

4=  Artillery  units (155mm) 

5=  Mechanized  battallion  (M-l’s) 

6=  Aircraft  Fighter  Squadron  (F-16’s) 

7=  Medium  truck  companies 
8=  USAF  Airlift  Control  Elements  (ALCE) 
Parameters  Used: 

ATTRIT(i,j):  the  attrition  rate  for  type  "i"  aircraft 
flying  on  mission  type  "j". 

ADC(i):  airdrop  capability  of  aircraft  type  "i"  in 
number  of  standard  pallets. 

AT(y):  anti-tank  capability  of  type  "y"  unit 
expressed  in  equivalent  TOEs. 

CI(y):  the  combat  indicator:  1=  combat  unit, 

-1=  combat  support,  0=  neither  combat  or  combat 
support  unit. 

CPI(l):  the  combat  power  index  for  arriving  at  the 
front  in  period  1. 

CPLOST(i):  the  expected  combat  power  lost  for  every 
aircraft  load  of  type  “ i"  lost  due  to  combat 


DAPDFT:  the  distance  between  the  APOD  and  the  FOL. 

DUSAPD:  the  distance  between  the  U.S.  and  the  APOD. 

DUSFRT:  the  distance  between  the  U.S.  and  the  FOL. 

EAS(i,k):  the  standardizing  factor  for  the  material 
handling  equipment  required  to  unload  aircraft 
type  " i“  with  type  ‘'k1'  cargo. 

FLT(y):  the  front  line  trace  capability  of  a  type 
"y"  unit. 

FP(y):  the  firepower  capability  of  type  "y"  unit. 

GAT(l):  the  desired  goal  for  anti-tank  power  by 
period  "1". 

GFLT(l):  the  desired  goal  for  front  line  trace  by 
period  "1". 

GCF(l):  the  desired  goal  for  combat  power  by  period 

"l". 

GT(i):  the  authorized  ground  time  for  type  "i" 
aircraft  to  offload/onload  and  refuel  if 
required  (hours/aircraft). 

Ha:  the  material  handling  equipment  capacity  in 

pallet  equivalents  existing  at  the  APOD  prior 
to  the  deployment  (pallets/day ) . 

Hf :  the  material  handling  equipment  capacity  in 

pallet  equivalents  existing  at  the  FOL  prior 
to  the  deployment  (pallets/day). 

INTER(i):  the  number  of  days  required  for  aircraft 
type  "i"  between  intertheater  missions  (days). 

INTRA(i):  the  number  of  days  required  for  aircraft 


,t*  ,t|  ,»»  ♦»  .»»  _♦**  **'  4+~'  tt’  It  *»' 


type  "i"  between  intratheater  missions  (days). 

Cargo(i,k):  the  cargo  capacity  of  type  “ i ”  aircraft 
with  type  "k“  cargo  (tons/ac). 

MHE(i):  the  pallet  capacity  of  type  "i"  aircraft  for 
airland  missions  (463L  pallets). 

My:  the  movement  time  required  by  a  unit  to  go  from 
the  APOD  to  the  FOL  by  ground  transportation, 
rounded  to  the  nearest  integer  (periods). 

L:  the  number  of  periods  considered  in  the  model. 

For  this  model  four  periods  are  considered,  each 
period  lasting  five  days. 

NPAL:  the  number  of  pallet  equivalents  that  an  ALCE 
unit  can  manage  in  one  day  (pallets/day) . 

NPRKA(i):  the  number  of  type  ”i“  aircraft  that  can 

simultaneously  park  on  the  APOD  ramp  (aircraft). 

NPRKF(i):  the  number  of  type  " i"  aircraft  that  can 

simultaneously  park  on  the  FOL  ramp  (aircraft). 

NPK(y):  the  number  of  type  "y"  unit  fighter’s  that 
can  simultaneously  park  on  the  APOD  ramp  (ac). 

NTAC(y):  the  number  of  fighter  aircraft  assigned  to 
type  "y"  unit  (aircraft/unit). 

NUM(i,l):  the  number  of  type  ”i"  aircraft  available 
for  the  first  period  and  committed  for  this 
specific  operation  (aircraft). 

NUM(i,l):  the  number  of  aircraft  available  for  period 
”1" (after  attrition  is  considered). 

PL:  the  number  of  days  in  one  period.  For  this  model 
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PL  =  5  days. 


RC:  the  number  of  pallets  the  Army  Riggers  can 
prepare  in  one  day  (pallets/day ) . 

SPD(i):  the  no-wind  cruise  speed  of  type  "i"  aircraft 

TON ( y , k ) :  the  amount  of  tonnage  that  must  be  moved 
for  type  "y"  unit  of  type  "k"  cargo  (tons/ac). 

T0N(y,3):  the  supply  consumption  rate  of  type  "y" 
unit  (tons/unit/day ) . 

TMC(y):  the  ton-mileage  capability  for  the  transpor¬ 
tation  of  supplies  to  the  front  for  type  "y" 
units  (ton/miles/unit/day ) . 

TVL(y):  the  distance  that  unit  type  "y"  can  travel 
in  one  day. 

Onit(y):  the  total  number  of  type  "y"  units  available 

UTE(i,j):  the  utilization  rate  of  type  “i"  aircraft 
on  a  type  "j”  mission  (hours/day). 


ObJostivo 

When  attempting  to  maximize  the  amount  of  combat 
power  delivered  to  a  theater  commander,  it  is  evident  that 
the  scenario  will  affect  which  types  of  units  will  be  more 
effective.  The  three  criteria  used  in  this  model  to 
quantify  combat  power  are  anti-tank  capability,  firepower, 
and  defense  frontage  (front  line  trace  or  FLT)  (7:52). 

All  three  measures  of  effectiveness  contribute  to  total 
combat  power  delivered,  but  each  would  possess  a  certain 


advantage  given  different  scenario  limitations. 

Otherwise,  one  could  state  that  the  combat  attributes  were 
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"mutually  reinforcing"  and  attempt  to  maximize  all  three 
variables  regardless  of  the  types  of  units  deployed  or 
scenario  parameters.  In  reality,  if  the  scenario  involves 
a  large  mechanized  enemy  force,  the  deployed  units’ 
ability  to  project  a  strong  anti-tank  capability  may  well 
override  its  ability  to  defend  a  given  amount  of  ground. 
These  real-world  constraints  dictate  a  model  which  allows 
the  user  the  opportunity  to  stress  the  importance  of  one 
measure  of  combat  power  over  another. 

Two  of  these  three  measures  of  combat  power  are 
independent  of  time,  whereas  the  third  is  very  time 
dependent.  The  amount  of  anti-tank  capability  and 
defensive  frontage  are  independent  of  time.  The  third 
measure  of  combat  power,  firepower,  emphasizes  a  point 
discussed  in  the  unclassified  portion  of  the 
Congress ionally  Mandated  Mobility  Study  (8:30).  This 
concept  stressed  the  importance  of  timely  delivery  and 
stated  that  the  delivery  of  a  small  force  to  the  correct 
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location  can  be  up  to  six  times  as  effective  as  the  same 
force  delivered  late.  The  time  independent  goals  of  anti¬ 
tank  capability  (AT)  and  defensive  frontage  (FLT)  will  be 
discussed  first.  Anti-tank  capability  is  represented  by 
using  the  multiplier  AT(y)  for  each  of  "y"  units.  For 
time  period  1’,  total  anti-tank  potential  is  described 
with  the  following  equation: 
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Y’  3  1J 

l  l  Z  AT(y)*u(y,m, 1) 

y-1  ra=l  1=1 

This  general  expression  must  now  be  broken  down  into 
the  different  methods  used  to  deliver  units  to  the  front. 
As  discussed  in  an  earlier  section  of  this  chapter,  m=l  i 
used  to  represent  direct  delivery  to  the  front;  m=2 
represents  self-propelled  units  which  arrive  at  the  APOD 
then  travel  to  the  forward  operating  location  on  their 
own.  The  capability  of  these  units  should  not  be 
considered  until  the  units  actually  arrive  at  the  front 
and  the  units’  assets  are  available  to  the  area 
commander.  To  consider  the  delay,  the  variable  "My"  is 
defined  as  the  number  of  periods  required  by  the  unit  to 
travel  from  APOD  to  FOL.  The  constraint  for  m=2  now 
becomes : 

Y’  l’-My 

E  2  AT ( y )  *  u(y , 2 , 1 ) 
y=l  1=1 

The  representation  for  units  requiring  intratheater 
airlift  is  (m=3).  As  stated  by  Cooke,  the  time  required 
to  go  from  APOD  to  FOL  is  not  considered  because  the 
enroute  time  is  signif icantly  less  than  the  time  for  one 
period  to  elapse.  The  expression  becomes: 

Y’  1’ 

2  S  AT ( y )  *  u(y , 3 , 1 ) 
y= 1  1=1 
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The  last  area  to  be  considered  consists  of  units 
assigned  to  the  APOD  that  possess  an  anti-tank  potential. 
The  next  equation  considers  units  at  the  APOD  that 
contribute  to  the  total  anti-tank  power: 

J/TC7)  *  u(jr,2.1) 

Summing  the  possible  sources  of  anti-tank  power  for 
period  1*  produces  the  following  equation: 

Y  1*  Y’  l'-My 

Si  EiAT<y)  *  u(y,m,l)  +  2  2  AT(y)  *  u(y,2,l)  + 

y=l  1=1  m=l , 3  y- 1  1=1 

Y  1’ 

S  ,  S  AT(y)  *  u(y,2,l)  2  GAT(l’)  (1) 

y=Y  +1  1=1 

The  first  summation  encompasses  direct  (m=l)  and 
intratheater  delivery  (m=3),  while  the  second  term  (m=2) 
includes  units  that  reach  the  front  via  ground 
transportation,  and  the  third  summation  includes  units 
assigned  to  the  APOD.  ( 7 : 53 ; 33 : Chap  3,8). 

Defensive  frontage  (FLT)  is  developed  in  a  similar 
manner.  In  this  case,  the  multiplier  FLT(y)  is  used  to 
represent  the  amount  of  defensive  frontage  capable  of 
being  defended  by  a  specific  unit.  The  development  of  the 
constraint  follows  the  same  format  as  the  previous 
constraint  and  results  in  the  following  equation  for 
period  1 ’ : 


nnHmaiT 


III  VI  J1 


APOD.  The  second  summation  includes  the  requirement  to 
deliver  the  men  and  equipment  to  the  objective  area  before 
the  combat  power  index  is  applicable.  The  third  summation 
includes  the  firepower  deployed  to  the  APOD  which  is  not 
to  be  transported  to  the  front.  ( 7 : 55 ; 33 : Chap . 3 , 9-10 ) . 

Any  of  these  three  measures  of  effectiveness  (MOE’s) 
can  be  used  interchangeably  as  the  objective  function 
based  upon  the  requirements  of  the  scenario.  The  other 
two  equations  would  then  be  entered  in  the  constraint 
matrix.  These  three  equations  represent  the  combat  power 
capability  of  the  units  to  be  airlifted  during  the 
contingency.  Only  the  equation  on  firepower  delivered  is 
t ime- dependent ,  and  that  concept  is  demonstrated  by  using 
the  CPI(l)  time  delay  The  first  two  equations  are  not 
dependent  upon  time,  and  measure  combat  power  based  upon 
the  unit  type  and  the  quantity  of  that  unit  delivered  to 
the  objective  area. 

An  examination  of  the  three  equations  reveals  a 
consistent  relative  ranking  of  the  units.  Regardless  of 
whether  anti-tank  power,  defensive  frontage,  or  fire  power 
is  maximized,  the  same  units  will  be  delivered  to  the 
objective  area.  Units  of  type  y=4  possess  the  largest 
percentage  of  fire  power,  defensive  frontage,  and  anti¬ 
tank  power.  The  relative  relationship  between  the  units 
is  preserved  with  each  of  the  three  possible  goals.  By 
simplifying  the  scenario  so  the  objective  is  to  maximize 
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fire  power,  the  other  two  measures  of  effectiveness  can  be 
set  as  constraints.  As  long  as  these  two  constraints  are 
set  at  a  feasible  level,  the  objective  function  of 
maximizing  fire  power  is  left  relatively  unaffected  by  the 
setting  of  defensive  frontage  or  anti-tank  capability 
levels. 

For  the  purposes  of  this  scenario,  the  goal  is  to 
maximize  fire  power  delivered  to  the  theater  commander, 
and  the  two  other  measures  of  effectiveness  can  be  set  as 
greater  than  or  equal  to  constraints.  The  objective 
function  equation  becomes: 

MAXIMIZE: 

Y*  1»  . 

E  2  Z  CPI(l)  *  FP(y )  *  u(y ,m, 1 )  + 
y=l  1=1  m=l,3 

Y*  i’-Mx 

E  E  CrI ( 1  +  My)  *  FP(y)  *  u(y,2,l)  + 
y=l  1=1 

Y  1*  I 

S  Vf  ,  E  CPI(l)*FP(y)*u(y,2,l)-E  CPLOSTC i )*ATT( i )  (4) 

y=Y  +1  1=1  i=l 

The  next  section  of  this  chapter  discusses  the 
aircraft,  cargo,  and  airfield  restrictions  that  constrain 
the  maximization  of  combat  power  delivered  to  the  field 
commanders.  This  section  also  further  defines  the  terms 
CPLOST(i)  and  ATT(i)  in  the  equation  above. 


The  quantity  of  combat  power  delivered  to  an 
objective  area  by  air  is  directly  affected  by  the 


characteristics  of  the  units  to  be  deployed  and  the 
capabilities  of  the  aircraft  involved  in  the  airlift. 

Each  of  the  constraints  to  follow  will  be  discussed  in 
detail  in  order  to  examine  the  influence  that  specific 
constraint  has  upon  the  maximization  of  combat  power 
delivered  to  the  theater  commanders  over  a  given  period  of 
time.  The  equations  were  formulated  largely  by  Cooke, 
updated  by  Tate,  and  changes  were  made  here  when  necessary 
to  reflect  conceptual  iifferences  or  changes  in  the  data 
base. 

Aircraft Limitations 

The  time  required  to  move  units  from  the  CONUS  to  the 
objective  area  is  heavily  influenced  by  the  total  number 
of  aircraft  dedicated  to  the  deployment  and  the  ability  of 
each  aircraft  type  to  fly  at  a  sustained  flying  hour 
level.  The  first  constraint  deals  with  the  number  of 
aircraft  sorties  as  determined  by  the  number  of  aircraft 
available  and  the  time  required  for  an  aircraft  to  fly  a 
complete  mission.  A  complete  mission  involves  departing 
home  station,  flying  the  required  mission  legs,  and 
returning  to  home  station  after  standard  ground  times,  and 
expected  maintenance  down  times,  ready  for  another  home 


station  departure.  The  second  constraint  determines  the 
expected  number  of  aircraft  sorties  capable  of  being 
generated  by  each  aircraft  type  (DTE  rate),  based  upon  the 
number  of  hours  each  aircraft  can  be  expected  to  fly. 

Aircraft  Sorties.  The  first  constraint  determines 
the  number  of  aircraft  sorties,  by  aircraft  type,  based 
upon  the  number  of  aircraft  expected  to  be  available  for 
the  deployment  of  the  combat  units.  The  total  number  of 
sorties  is  also  dependent  upon  the  time  required  for  an 
aircraft  to  depart  home  station,  fly  the  proposed  routing 
with  standard  ground  times  and  expected  maintenance  down 
times,  and  return  to  home  station  ready  for  smother 
mission.  The  range  and  speed  of  each  aircraft  differ 
slightly,  and  for  this  specific  scenario,  data  was 
obtained  from  MACR55-141  and  HQ  MAC/XOSS  (12).  HQ  MAC 
provided  the  estimated  distances  to  the  various  bases  and 
forecasted  wind  factors  for  the  proposed  routing. 

Estimated  ground  speeds  of  the  various  aircraft  were 
determined  based  upon  true  airspeeds  and  wind  factors  for 
each  mission  leg.  Using  the  estimated  distances  and 
groundspeeds ,  an  enroute  time  was  determined  for  each  leg 


of  the  mission.  These  estimated  enroute  times  were  added 
together  and  combined  with  the  estimated  maintenance 
downtimes  and  scheduled  ground  times  for  each  aircraft 
type  to  determine  a  total  time  required  to  be  ready  for 
another  mission. 
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Cooke  assumed  overflow  between  days  because  his 
period  was  equal  to  one  day,  yet  the  mathematical 
formulation  used  period  lengths  of  five  days.  Using  a 
similar  time  frame  of  four  linked  five-day  periods  will 
adequately  represent  the  time  period  considered  most 
critical  for  testing  the  effectiveness  of  airlift  in 
delivering  combat  power  to  a  given  objective  area. 

Thus  far  two  factors  have  been  discussed  which  affect  the 
total  number  of  aircraft  sorties  available  for  the 
airlift,  the  time  required  for  an  aircraft  to  transit  the 
proposed  route  structure  and  the  number  of  aircraft 
expected  to  be  available  for  the  operation.  A  third 
factor  affecting  sortie  generation  is  possible  aircraft 
attrition.  Attrition  would  be  related  to  the  aircraft  and 
mission  type.  Attrition  for  a  specific  aircraft  and 
mission  type  is  a  subjective  determination  with  a  number 
of  possible  complicating  factors.  While  some  theorists 
believe  a  direct  delivery  aircraft,  such  as  the  C-17,  is 
more  susceptible  to  enemy  fire  because  it  flies  closer  to 
the  actual  battle  area,  other  sources  believe  that  direct 
delivery  disperses  the  airlift  force  over  a  wider  area, 
thereby,  enhancing  the  survivability  of  all  airlift 
aircraft  operating  in  a  threat  environment.  (30:33).  The 
C-17  will  possess  increased  capability  to  airland  and 
airdrop  men  and  supplies  directly  to  the  objective  area 
while  minimizing  the  threat  to  the  aircraft,  but  the 
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threat  still  exists.  A  further  discussion  of  attrition 


follows  in  the  development  of  the  scenario.  Attrition 
variables  are  included  in  the  model  and  the  coefficients 
have  been  determined  subjectively.  The  purpose  of 
including  attrition  in  the  model  is  not  to  validate  the 
attrition  values  but  to  examine  the  effects  of  attrition 
on  the  delivery  of  combat  power.  Attrition  values  will  no 
doubt  be  debatable,  but  hopefully  the  relative  scaling  of 
the  attrition  values  are  reasonable  estimates  of  potential 
attrition  in  a  low  threat  environment.  The  mathematical 
equation  depicts  the  negative  effect  of  attrition  on 
sorties  available.  The  attrition  of  the  number  of 
aircraft  available  for  a  given  period  of  time  is  developed 
by  first  setting  the  expected  number  of  aircraft  available 
in  the  first  period  with  the  variable  NUM(i,l).  Aircraft 
type  "i"  begins  the  operation  in  subsequent  periods  with  a 
given  number  of  aircraft  available  from  period  one  minus 
any  attrition  of  aircraft.  The  relationship  is  NUM(i,l)  = 
NUM(i,l)  -  Attrition.  Attrition  depends  upon  aircraft 
type  and  mission  type  and  is  formulated  as: 

r-i  j  k 

Attrition  =  S  S  S  Attrit(i,j)  *' x(i , j , k, 1 ’ ) 

1=1  j=l  k=l 

Since  we  are  dealing  with  periods,  and  not  days  the 
formulation  becomes: 


For  1=1:  NUM(i,l)*Pl  =  NUM(i,l)*PL 


w  imiwvw*  iMijmsscmwwr 


rxiiunim 


nnrnniimnrvimKramanwmamjrTaiw  wwtoj 


For  1=2, 3, 4: 


NUM( i , 1 ) *P1=NUM( i , 1 ) *P1 


l’-l  J  K 

-2  2  2  Attrit ( i , j ) *PL  *  x(i,j,k,l’) 

1=1  j=l  k=l 


By  rearranging  the  terms,  the  final  equation  to  determine 
the  number  of  sorties  available  during  the  deployment 
phase  is: 


For  1—  1  •  g  ^ 

Inter ( i )*  2  2  x(i,j,k,l)  +  Intra( i ) *  2  x( i , 4 , k, 1 ) 

j=l  k=l  k= 1 


<.  Num(  i  ,  1 )  *PL 


For  1=  2,3,4: 


3  K  K 

Inter(i)*  2  2  x(i,j,k,l)  +  Intra(i)*  2  x(i,4,k,l)  + 

j=l  k=l  k=l 


l’-l  J  K 

2  2  2  Attrit ( i ,j)*x(i,j,k,l’ )  iNum(i,l)*PL  (5) 

1=1  j=l  k= 1 


There  will  be  one  constraint  for  each  aircraft  type  per 
period,  or  a  total  of  (I*L)  equations.  (33: Chap  3,11) 


(JTE  Rate.  The  second  set  of  constraints  related  to 
aircraft  limitations  consider  the  expected  usage  rate  of 
the  aircraft.  The  expected  utilization  (UTE)  rate  is  a 
commonly  used  method  of  measuring  the  airlift  capability 
of  an  aircraft.  UTE  rate  is  computed  by  dividing  the 


total  programmed  flying  hours  for  a  given  aircraft  type 
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(TPFH)  by  the  primary  aircraft  authorized  (PAA)  and 
dividing  that  number  by  360  to  get  an  expected  daily  usage 
rate  (12).  An  approximation  of  the  expected  surge 
capability  UTE  rate,  which  is  higher  than  the  standard 
"how  much  are  we  flying  today  rate"  is  used  for  this 
scenario.  The  average  number  of  sorties  each  particular 
aircraft  type  can  be  expected  to  fly  is  depicted  as 
follows : 

[PL*UTE(  i  ,  j  ) *Num(  i  ) *SPD(  i  )  ] /( 2*DUSAPD)  =  aircraft/period 

The  period  length  *  the  UTE  rate  *  number  of  aircraft 
available  *  speed  of  the  aircraft  divided  by  two  *  the 
distance  from  the  United  States  to  the  APOD  =  required 
aircraft/period.  As  noted  by  Cooke  and  Tate,  the  two  in 
the  denominator  is  used  to  denote  the  round  trip  distance 
involved  in  the  airlift.  The  terms  can  be  rearranged  to 
develop  the  following  constraint: 

J  K 

.2,  S  Cl/UTE(i, j)]*x(i, j,k,l' )  1  [PL*Num( i )*SPD( i ) ] / 
j=l  k=l 

( 2*DUSAPD) 

A  standardizing  factor  was  developed  by  Cooke  to 
convert  aircraft  of  the  same  type,  flying  different 
missions  to  an  equivalent  number  of  sorties  for  a 
respective  mission.  The  ratio  [  (distance  from 
other)/(distance  from  US  to  APOD)]  shows  that  there  are 
more  sorties  available  for  intratheater  (j=4)  missions 
than  for  intertheater  missions  (j=l,2,3).  It  should  also 


be  noted  that  any  aircraft  used  strictly  for  intratheater 
missions,  such  as  the  C-130,  should  be  considered 
separately.  The  number  of  sorties  they  can  generate  are 
dependent  upon  the  number  of  aircraft  available  and  the 
distance  from  the  APOD  to  the  FOL. 


The  general  expression  for  the  constraint  for  period 
1’  and  for  aircraft  types  1  through  5  is  as  follows: 


K 

2 

k=l 


2  l/UTE(i,j)* 

j  —  2  i  o 


DUSFRT 

DUSAPD 


K 

*x(i, j,k,l’ )  +  2  l/UTE(i,l)* 

k=l 


K  DAPDFT 

x(i, l,k, 1* )  +  2  1 /DTE i  i ,  4 )  *  -------  *x(i,4,k,l’)  i 

k=l  DUSAPD 

PL*Num( i )*SPD( i ) /( 2*DUSAPD)  (6) 


The  first  summation  accounts  for  the  direct  (j  =  2)  and 
airborne  (j=3)  missions,  the  second  summation  accounts  for 
intertheater  missions,  and  the  third  summation  accounts 
for  aircraft  capable  of  also  flying  intratheater  missions 
(j=4).  The  C-130  sorties,  used  strictly  for  intratheater 
missions,  use  the  following  equation: 

K 

2  1 /UTE (4,4)*x(6,4,k,l’ )  i  PL*Num(6)*SPD(i)/2*DAPDFT  (7) 
k=l 


There  is  one  constraint  for  each  type  of  aircraft  per 


period  for  a  total  of  (I*L)  constraints  (7:59-60;33: 
Chap3, 11-13) . 


Airport  facilities  provide  three  additional 
constraints  on  strategic  mobility.  The  first  constraint 
is  on  the  capacity  of  the  APOD  ramp  to  park  available 
aircraft,  the  second  is  the  capacity  of  the  FOL  ramp  to 
park  available  aircraft,  and  the  third  constraint  is  the 
capability  of  the  material  handling  equipment  (MHE)  to 
offload  arriving  aircraft  and  onload  aircraft  for 
departure.  As  in  previous  thesis  efforts  on  this  subject, 
a  linear  relationship  is  assumed  between  the  area  required 
to  park  different  aircraft  types.  As  an  example,  if  the 
ramp  space  is  sufficient  to  park  30  C-130’s  or  20  C-141’s, 
it  is  also  Sufficient  to  park  15  C-130’s  and  10  C-141’s  or 
any  other  convex  combination.  As  pointed  out  by  Cooke, 
such  approximations  are  sufficient  for  all  but  the 
smallest  of  airfields,  and  small  airfields  are  not  usually 
selected  as  main  operating  bases  because  of  ramp  space 
limitation  (7:61). 

Parking  space  per  unit  time  for  each  aircraft  type  is 
used  as  the  parameter  to  measure  available  ramp  space. 
Using  the  scheduled  ground  time  [GT(i)]  and  the  number  of 
aircraft  of  a  given  type  that  can  park  at  the  APOD 
[NPRK(i)],  it  is  possible  to  determine  how  many  sorties  of 
a  given  aircraft  type  can  be  generated  through  the  APOD 
in  a  single  period.  The  parking  space  is  constrained  as 


follows : 


I  K 

.2  2  ,  £  x( i , j ,k, 1)  *  GT( i )  /  [PL*NPRKA( i )*24]  <  1 

i=l  j  =  l ,  4  k=l 

When  the  left  hand  side  of  the  constraint  equals  one, 
ramp  space  on  the  field  is  completely  saturated.  The 
number  of  aircraft  for  each  type  that  can  fit  on  the  ramp 
was  obtained  from  a  HQ  MAC  maximum  on  the  ground  listing 
(MOG).  The  listing  includes  several  hundred  bases 
currently  surveyed  by  MAC  and  published  by  HQ  MAC/DEP. 

The  possibility  exists  that  fighter  aircraft  may  be 
assigned  to  the  APOD.  If  that  were  the  case,  they  would 
also  require  parking  spots  and  the  constraint  would  now 
be : 


K 


1 

2  -S,  .  2  x(i,j,k,l)  *  GT(i)/[PL*NPRKA(i)*24]  + 

i=l  j=l , 4  k=l 


£  ,  0 ( y , 2 , 1 ) *NTAC ( y ) /PL*NPK ( y )  £  1  (8) 

y=y  +1 


The  first  summation  includes  aircraft  flying  directly 
to  the  APOD  or  flying  intratheater  missions  between  the 
APOD  and  FOL.  The  second  summation  includes  any  fighter 
units  stationed  at  the  APOD.  There  is  one  constraint  per 
period  for  a  total  of  (L)  constraints  (7:62-63). 
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The  second  constraint  is  very  similar  to  the  first, 
except  it  determines  the  parking  capacity  at  the  FOL.  It 
was  not  included  in  previous  efforts  but  should  be  added 
for  completeness .  The  mathematical  formulation  is  as 
follows : 

4  K 

2  o  o  *  S  x(i, j,k,l)*GT(i)/[PL*NPRKF(i)*24]  i  1  (9) 

i  =  2 , 3 , 6  j-2  k=l 

The  equation  only  considers  C-17,  C-141  and  C-130’s 
because  the  other  aircraft  types  will  not  transit  the 
front.  As  in  the  previous  equation,  there  should  be  (L) 
constraints . 

The  third  constraint  determines  the  airfield’s 
capability  to  offload  and  onload  equipment  at  the  APOD  and 
offload  equipment  at  the  front.  The  time  required  to 
upload  or  download  an  aircraft  is  dependent  upon  the 
number  of  pallets  an  aircraft  can  carry  and  the  type  of 
cargo  to  be  moved.  The  number  of  pallets  (NPAL)  that  the 
various  aircraft  in  this  scenario  can  hold  is  listed  in 
Table  2.  The  ease  with  which  each  aircraft  can  be  loaded 
or  unloaded  [EAS(i,k)]  is  a  subjective  determination  based 
upon  the  aircraft  type,  cargo  type,  and  sophistication  of 
loading  equipment  (k-loaders  vs  fork  lifts).  If  all  the 
aircraft  can  be  loaded  with  the  same  degree  of  difficulty, 
all  aircraft  will  be  given  an  index  of  one.  The  arriving 
aircraft  can  expect  to  be  downloaded  by  the  equipment  (K- 
loaders,  forklifts,  etc.)  that  are  currently  available  at 


t. 


the  APOD  (Ha)  and  the  front  (Hf).  Once  an  ALCE  unit  is  in 
place,  it  will  provide  additional  offloading  capability. 

A  multiplier  is  used  to  approximate  the  additional 
assistance  gained  from  any  arriving  ALCE  unit.  For  this 
model  a  multiplier  of  0.5  is  selected  as  a  representative 
level  of  improvement  expected  by  the  arrival  of  ALCE 
units.  The  constraint  must  also  consider  the  restriction 
that  the  number  of  aircraft  moving  through  the  APOD  or  FOL 
cannot  be  greater  than  the  capabilities  of  the  airfield 
MHE  to  upload  or  download  them.  The  constraints  are  as 
follows : 


APOD: 

1  K 

2  E  S  EAS(i,k)*MHE(i)  *  x(i,j,k,l’) 

1=1  j=l , 4  k=l 

l’-l 

-  Z  NPAL*PL*u(8,2,l)-0.5*NFAL*PL*u(8,2,l’ )  <  Ha*PL  (10) 


E 

i=2 ,3,6 


E 

j  =  2,4 


K 

£  EAS ( i , k )  *  MHE(i)  *  x(i,j,k,l’) 
k- 1 


1’-  My 

-  2  E  NPAL*PL*u  (  8  ,  m ,  1 )  Hf  *PL  (11) 

1  =  1  oi=  1 1  3 


There  is  one  constraint  per  period  for  both  the  APOD 
and  the  FOL,  or  (2*L)  constraints,  ( 7 : 61-64 ; 33 : Chap  3,  13- 


This  constraint  encompasses  the  restriction  that  the 
number  of  units  shipped  must  be  equal  to  or  less  than  the 


total  number  of  units  available.  The  constraint  is 
expressed  as : 


3  L 

2  S  u(jr,m,l)  A  Unit ( y )  (12) 

m=l  1=1 

There  is  one  constraint  for  each  type  of  unit,  for  a 
total  of  Y  constraints.  ( 7 : 64 ; 33 : Chap  3,17) 

Shipment _af  Units 

Units  can  be  delivered  to  the  APOD  using  three 
primary  methods,  1)  via  the  APOD,  2)  via  direct  delivery, 
and  3)  aerial  delivery. 

The  first  section  discusses  the  delivery  of  units  to 
the  APOD.  When  discussing  the  delivery  of  units  to  a 
given  area,  the  entire  unit  must  be  delivered  before  the 
unit  is  considered  "closed".  Each  unit  must  be  delivered 
in  its  entirety.  Outsized  cargo,  oversized  cargo,  bulk 
cargo,  and  personnel  must  all  be  delivered  before  the  unit 
is  considered  ready  for  its  combat  or  support  mission. 
Excess  cargo  [A(k,l)]  or  early  cargo  [A(K,1-1)]  also  needs 


to  be  considered  in  formulating  the  constraint.  The 
complete  formulation  is  as  follows: 


1  Y  3 

2  cargo(i,k)  *x(i,l,k,l)-2  2  TON(y,k)  *  u(y,m,l) 

i  =  l  y=l  m=2 

+  A(k, 1-1 )  -  A(k, 1 )  =  0  (13) 


To  insure  that  cargo  moved  to  the  APOD  is  delivered 
to  the  front,  the  following  constraint  is  added: 


1  Y’ 

2  cargo(i.k)  *  x(i,4,k,l)  -  2  TON(y , k)*u(y , 3 , 1 ) 

i=l  y- 1 

+  P( k , 1-1 )  -  P(k, 1 )  =  0  (14) 

There  are  two  constraints  for  every  cargo  type  and  period 
for  a  total  of  (K*L)  constraints .( 7 : 65-66 ; 33 : chap  3,17-18) 

The  second  section  discusses  the  delivery  of  units  to 
the  front  via  direct  delivery.  The  mathematical 
expression  is  very  similar  to  the  APOD  equation  with  the 
exception  that  only  aircraft  flying  direct  missions  will 
be  considered: 


1  Y 

2  cargo(i.k)  *  x(i,2,k,l)  -  2  TON(y,k)*u(y , 1 , 1) 

i=l  y=l 

+  P(k, 1-1 )  -  P( k, 1 )  =  0  (15) 
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Again,  there  is  one  constraint  for  each  cargo  type  and 
period  for  a  total  of  (K*L)  constraints.  ( 7 : 66 ; 33 : chap  3 
18) 


The  last  method  to  be  discussed  is  the  use  of 
airdrop  missions  to  deliver  the  men  and  supplies.  An 
indicator  variable,  AB(y)  =  1  or  0,  is  used  to  indicate 
whether  the  unit  is  airdrop  capable  or  not.  The  only 
aircraft  considered  are  those  aircraft  capable  of 
intertheater  airdrop  missions,  the  C-17  and  the  C-141. 
While  the  C-130  is  also  airdrop  capable,  it  is  normally 
used  for  intratheater  missions  (no  air  refueling 
capability),  and  intratheater  airdrop  missions  are  not 
considered  in  this  scenario.  The  expression  is: 

Y’ 

E  cargo ( i,k)*x(i,3,k,l)  -  E  AB( y ) *T0N(y , k) *u( y , 1,1) 
i  =  2 , 3  y=l 

+  P(k, 1-1 )  -  P(k, 1 )  =  0  (16) 

As  stated  before,  there  is  one  constraint  for  each 
period  and  cargo  type,  or  (L*K)  constraints.  ( 7 : 66 ; 33 : Chap 
3,  18-19). 

Shipmenti-of  Supplies 

Supplies  are  considered  bulk  items  (K=3),  and  once  a 
unit  is  delivered  to  the  objective  location,  supplies  must 
be  kept  at  a  stabilized  level  in  order  to  sustain  the 
unit’s  operation.  The  input  of  supplies  must  exceed  the 


consumption  of  supplies  for  all  the  deployed  units.  The 
shipment  of  supplies  will  be  discussed  in  three  separate 
sections.  The  first  section  discusses  the  requirement 
that  supplies  reach  the  units  at  the  front,  and  the  last 
section  discusses  the  need  to  adequately  supply  units 
based  at  the  APOD. 

This  section  discusses  the  need  to  insure  adequate 
supplies  reach  the  front.  An  assumption  for  this 
constraint  is  that  any  unit  not  based  at  the  APOD  is  based 
at  the  FOL  and  will  require  its  supplies  to  reach  the 
front.  Units  such  as  the  fighter  squadron  and 
transportation  units  are  expected  to  require  APOD 
supplies;  whereas,  the  other  units  in  the  model  will 
require  supplies  at  the  FOL.  It  is  assumed  that  units 
arriving  in  a  given  period  have  enough  supplies  to  sustain 
the  unit  for  the  remainder  of  the  period.  With  this 
assumption,  the  requirement  is  to  insure  adequate  supplies 
are  delivered  to  units  arriving  in  previous  periods  only. 
As  stated  in  Tate’s  thesis,  supplies  can  be  delivered  to 
the  front  via  direct  delivery,  airborne  delivery, 
transportation  via  intratheater  airlift,  and  overland 
shipment  by  truck.  Leaving  the  truck  shipments  to  be 
discussed  later,  the  total  amount  of  supplies  delivered  by 
airlift  is  found  by  adding  the  supplies  delivered  by 
direct  delivery,  (j=2),  airborne  delivery  (j=3),  and 
intratheater  transshipment  (j=4).  The  supplies 
transported  are  represented  by: 
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1  4 

2  . 2  cargo ( i , 3 )  *  x(i,J,3,l’) 
i=l  j  =  2 

For  those  units  being  shipped  by  truck,  it  is 
necessary  to  consider  the  time  required  to  travel  from  the 
APOD  to  the  FOL  and  the  capability  of  the  unit  to  move 
supplies.  The  parameter  TRP(y)  is  used  to  designate  the 
capability  of  a  unit  to  move  supplies.  This  parameter  is 
applicable  to  units  such  as  a  truck  company  which  can  be 
expected  to  deliver  supplies  from  the  APOD  to  the  FOL. 

The  constraint  expressing  the  number  of  tons  a  unit  can 
move  in  one  period  is: 

TRP ( y ) *PL/2*DAPDFT 

The  two  in  the  denominator  reflects  the  round  trip 
distance  from  the  APOD  to  the  front.  The  total  tonnage 
shipped  between  the  APOD  and  FOL  by  truck  is  as  follows: 

Y  3  1 ’ -My 

S4  E  E  TRP(y)*PL/2*DAPDFT]*u(y ,m, 1) 
y=l  m=l  1=1 

The  next  equation  combines  the  expressions  for 
delivering  supplies  to  the  front  and  subtracting  the 
consumption  of  the  units  previously  delivered  to  the 
front.  The  equation  is: 


14  Y  3  l’-My 

S  E  cargo(i,3)*x(i, j,3,l’ )  +  2  E  2  [TRP(y)*PL/ 
i=l  j  =  2  y= 1  m=l  1=1 


Y’  3  l’-l 

2*DAPDFT]*u(y , m, 1 )  -  2  2  2  PL*TON( y , 3 )*u(y , m , 1 )  2  0 

y=l  m=l  1=1 


An  additional  constraint  must  be  added  to  insure  that 


the  aircraft  and  trucks  do  not  attempt  to  deliver  supplies 


that  are  not  available  at  the  APOD.  Letting  SI  equal  the 


supplies  transshipped  during  period  1  and  formulating  the 


constraint  such  that  the  amount  of  supplies  delivered  by 


airlift  plus  the  supplies  delivered  by  truck  are  greater 


than  or  equal  to  the  supplies  delivered  during  the  period, 


the  expression  for  period  1’  becomes: 


I  Y  3  l’-My 

Ecargo(i,3)*x(i, 4,3,1’ )  +  2  2  2  [TRF(y)*PL]/ 

1=1  y=l  m=l  1=1 


( 2*DAPDFT)  *  u ( y , m , 1 )  >  SI 


The  final  equation,  letting  the  parameter  T0N(y,3) 


represent  the  daily  consumption  rate  for  unit  "y", 


PS(3,1’)  represent  any  surplus  supplies  shipped  at  the  end 


of  the  1 ’ th  period,  and  PS(3,1’-1)  represent  the  surplus 


from  previous  periods,  is-' 


14  Y  3  1  ’  -My 

2  2  cargo ( i , 3 )  *  x(i,j,3,l’)  +  2  2  2  [TRP(y )  *  PL/ 

i=l  j=l  y= 1  m= 1  1=1 


Y’  3  l’-l 

2*DAPDFT]  *  u(y,m,l)  -  2  2  2  PL*T0N(y,3)  *  u(y,m,l) 

y=l  m=l  1=1 


+  SI  -  PS ( 3 , 1 ’  - 1 )  +  PS ( 3 , 1 ) 


0 


(  17) 


There  are  two  constraints  per  period  or  a  total  of  ( 2 


*L)  constraints  ( 7 : 68-9 ; 33 : Chap . 3 , 22-23 ) . 

The  next  equation  deals  with  the  constraint  of 
supplies  for  the  APOD.  As  in  the  previous  cases,  the 
supplies  reaching  the  APOD  can  be  no  less  than  the 
requirements  of  the  deployed  units.  All  the  supplies 
delivered  to  the  APOD  are  either  consumed  at  the  APOD, 
delivered  by  intratheater  airlift  to  the  front,  or  shipped 
by  trucks  to  the  front.  The  supply  constraint  is  the  sum 
of  the  supplies  delivered  to  the  APOD  for  this  period  plus 
any  surplus  supplies  from  the  previous  period,  minus  the 
amount  of  supplies  either  consumed  at  the  APOD, 
transshipped  to  the  front,  and  any  excess  supplies  for  the 
period.  The  equation  includes  the  same  SI  variables  which 
equates  to  the  amount  of  supplies  delivered  from  the  APOD 
to  the  front  for  a  given  period.  This  scenario  assumes 
that  the  battalion  headquarters,  F-16  squadron,  and  the 
truck  company  receive  their  supplies  at  the  APOD.  The 
resulting  equation  is: 

I 

^E^cargo( i , 3 )  *  x(i, 1,3,1’)  +  A(3,l’-1)  -  SI  -  A(3,l’) 


Y  1’  -1 

S  PL  *  TON(y , 3 )  *  u(y,2,l)  =  0  (13) 

y=Y  +1  1=1 
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For  clarification,  the  first  Y’  units  are  deployed  to 


the  front  while  the  remainder  are  deployed  at  the  APOD. 
There  is  one  constraint  per  period  for  a  total  of  (L) 
constraints.  ( 7 : 70-71 ; 33 : Chap  3,23-24). 


Onit  Linkage 

This  section  deals  with  the  relationship  between 
distinct  but  related  unit  types.  Certain  units  require 
that  a  minimum  number  of  other  related  units  be  deployed 
if  they  are  to  be  deployed.  This  relationship  is  examined 
in  the  model  by  requiring  a  linkage  between  the  number  of 
combat  units  and  combat  support  units,  such  as  the  link 
between  a  battalion  and  battalion  headquarters.  These 
linkages  establish  certain  ceilings  where  only  a  given 
amount  of  one  unit  can  be  deployed  unless  another  related 
unit  increases  its  deployed  forces.  The  inverse 
relationship  exists  for  floors  of  certain  units  whereby 
particular  units  must  be  deployed  in  order  for  some 
distinct  but  related  unit  to  increase  its  resources. 

As  an  example  of  a  ceiling,  the  following  equation 
requires  a  minimum  number  of  combat  support  units  to  be 
deployed  to  backup  the  combat  units.  A  certain  ratio 


is  established  between  combat  and  combat- support  units. 

The  next  equation  uses  a  multiplier  [CI(y)]  to  distinguish 
between  combat,  combat  support,  and  other  units.  The 
multiplier  is  given  the  value  of  1,-1,  and  0  respectively. 
An  "other"  unit  may  include  an  ALCE  unit  which  is  neither 
a  combat  unit  or  combat  support  unit.  The  formulation  is 
as  follows: 

Y  3  L 

Z  Z  2  Cl (y )  *  u( y , m, 1 )  >  0  (19) 

y- 1  m=l  1=1 

This  equation  adds  one  constraint  to  the  problem. 
(7:71,  33: Chap  3,24-25). 

A  floor  is  another  typical  unit  linkage  which  requires 
deployment  of  a  minimum  of  a  certain  type  of  unit  given 
that  other  units  are  already  deployed.  This  relationship 
exists  between  combat  units  and  combat  support  units,  such 
as  a  battalion  headquarters.  One  method  of  linking  a  one 
to  five  ratio  between  HQ  units  to  battalion  units  is-' 

combat  unit  -  5*Hq  unit  <.0  (20) 

This  formulation  forces  a  HQ  unit  to  precede  the 
combat  battalion  in  the  initial  deployment,  which  is 
neither  desirable  nor  realistic.  The  following 
expression:  combat-3HQ  _<  2,  allows  up  to  two  combat  units 


to  deploy  before  a  HQ  unit  is  required  and  allows  a 
reasonable  number  of  combat  units  to  deploy  for  every  HQ 
unit  deployed.  This  equation  requires  an  additional 
constraint  for  each  linkage  between  combat  units  and  other 
types  of  units  in  the  scenario.  ( 7 : 72 ; 33 : Chap  3,25). 


This  constraint  examines  the  limitation  on  airdrop 
capability  imposed  by  the  ability  of  rigging  personnel  to 
properly  configure  the  units’  assets  into  pallets  capable 
of  being  airdropped.  Riggers  can  only  rig  so  many  pallets 
in  a  given  time  period,  and  this  expression  can  be 
expressed  mathematically  as: 

I 

£  ADC (i,3)*x(i,3,3,l)  -  P(3,l-1)  +  P(3,l)  =  RC*PL  (21) 
i  =  l 

The  variable  ADC(i,3)  refers  to  the  number  of  pallets 
an  aircraft  (i)  can  hold  in  the  airdrop  configuration,  and 
RC  refers  to  the  rigging  capability  of  the  army  riggers  in 
pallets  per  day.  The  P(3,l-1)  and  P(3,l)  represent  the 
previous  periods  surpluses  and  the  present  periods 
surpluses.  There  is  one  constraint  per  period  or  (L) 
constraints  ( 7 : 73 ; 33 : Chap  3,25-26). 


The  purpose  of  this  chapter  is  to  examine  the 


mathematical  equations  used  to  derive  the 
interrelationships  between  aircraft  capabilities,  cargo 
characteristics,  and  airfield  restrictions.  The  majority 
of  the  equations  were  derived  by  Cooke,  and  with  only 
minor  changes,  the  variable  names  were  determined  by 
Tate.  The  data  base  that  is  applied  to  these  mathematical 
equations  varies  from  previous  efforts  in  numerous  areas 
and  will  be  discussed  in  the  next  chapter.  This  chapter 
concludes  with  a  summary  of  the  objective  and  constraints 
as  follows: 


OBJECTIVE: 


Fire  Power  Capability 


MAXIMIZE: 

y>  l*  Y’  l’-My 

2  2  2  CPI ( 1 ) *FP (y)*u(y,m,l)  +  2,  2  CPI ( 1+My )*FP(y ) * 

y= 1  1=1  m=l , 3  y=l  1=1 


Y  1’ 

u(y , 2 , 1 )  +2  2  CPI ( 1 )*FP(y ) *u(y ,2,1) 

y=Y  +1  1=1 


1 

2  CPLOST(i)  *  ATT(i)  (4) 


SUBJECT  TO  THE  CONSTRAINTS: 


Anti-Tank  Capability 


Y*  1’  Y*  l’-My 

£  £  2,  AT(y)*u(y,m, 1)  +  2,  2  AT(y)*u(y, 2, 1)  + 

y=l  1=1  m=l , 3  y=l  1=1 

Y  1* 

2  2  AT(y)*u(y, 2, 1)  2  Anti-tank  goal 

y=Y’+l  1=1 


(1) 


Front  Line  Trace  Capability 


Y’  1’  Y’  l’-My 

r-1  1=1  m=l  fLT(y)*u(y>m'1)  +ylx  1S1FLT(y)*u(y,2fl)  + 


Y 

r 


1* 


2  FLT(y)*u(y , 2 , 1 )  2  Front  Line  Trace  Goal  .  (2) 
y=Y  +1  1=1 


Aircraft  Sortie  Generation  Capability  (One  for  each 
aircraft  type  and  period) 


3  K  K 

Inter(i)*  ^2  ^x(  i  ,  j  ,  k,  1 )  +  Intra(  i  i  ,  4  ,  k,  1 )  + 

l’-l  J  K 

2  2  2  Attrit ( i,j)*x(i,j,k,l’ )  i  Num(i,l)*PL  (5) 

1=1  j= 1  k= 1 
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fctfw  H  ^  C  bdw  II  CJ  >-4Ui  II 


Aircraft.  Utilization  (UTE)  Rate  (One  for  each  aircraft 
type  and  period) 


1 


2  l/UTE(i,j)* 
3  =  2,3 


DUSFRT 

DuiAPD 


K 

*x(i, j.k.l* )  +  E  1 /UTE ( i , 1 ) * 
k=l 


i ,  1 , k, 1 ’  ) 


K 

+  £  l/UTE(i,4)* 

k=l 


DAPDFT 

DUSAPD 


*x( i , 4 , k, 1 ’  ) 


<  PL*Num( i )*SPD( i ) /( 2*DUSAPD)  (6) 


tratheater  (UTE)  Rate  (C-130) 


l/UTE(4,4)*x(i,4,k,l')  ^  PL*Num(i)*SPD(i)/2*DAPDFT  (7) 


OP  Ramp  Space  (One  for  each  period) 


K 

,  2  2  x( i , j , k, 1 )  *  GT(i)/[PL*NFRKA(i)*24]  + 

1  j  =  l,4  k=l 

Y 

2  u(y , 2 , 1 ) *NTAC ( y ) /PL*NPK( y )  £  1  (3) 

y=y  +l 


FPL  Ramp  Space  (One  for  each  period) 


4  K 

£  £  x( i , j , k, 1 ) *GT( i ) /[PL*NPRKF ( i ) *24 ]  £  1  (9) 

j=2  k= 1 


(One  for  each 


APOD  Material  Handling  Equipment  (MHE) 
period) 


K 


I 

S  2  4  2  EAS ( i , k ) *MHE ( i )  *x(i,j,k,l’) 

i=l  j=l , 4  k=l 


l’-l 

•  SNPAL*PL*u( 8,2,1 )-0.5*NPAL*PL*u( 8,2,1* )  i  Ha*PL  (10) 


FPL  Material  Handling  Equipment  (MHE) 
period) 


(One  for  each 


,  „  2  .  E  2  EAS(i.k)  *  MHE(i)  *  x(i,j,k,l*) 

i=2, 3, 6  j=2 , 4  k=l 


1’-  My 

2  2  NPAL*PL*u ( 8 , m , 1 )  4  Hf  *  PL 

1=1  m= 1 , 3 


(11 


Unit  Limitations  (One  for  each  unit  type) 


3  L 

2  2  u ( y , m , 1 )  i  Unit(y)  (12 

m=l  1=1 


Shipment  of  Units  To  the  APOD  (Two  for  each  period  and 
cargo  size) 


I  Y  3 

.2  cargo (i,k)  *x(i,l,k,;)-2  2  TON(y,k)  *  u(y.m.l) 

i=l  y=l  m=2 

+  A(k, 1-1 )  -  A(k, 1 )  =  0  (13 
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To  insure  that  cargo  moved  to  the  APOD  is  delivered 


to  the  front,  the  following  constraint  is  added: 

I  Y' 

i  cargo(i,k)  *  x(i,4,k,l)  -  E  TON(y , k)*u(y , 3,1) 

1=1  y~ 1 

+  P(k, 1-1 )  -  P(k,l)  =  0  (14) 

S.hi PfO^nt  <?f — Units  via  Direct  Delivery  (One  for  each  period 
and  cargo  size) 


I 

£  cargo(i,k) 
i  =  l 


Y 

*  x( i , 2 , k, 1 )  -  E  TON (y,k)*u(y,l,l) 
y- 1 

+  P(k,  1-1 )  -  P(k,  1 )  =  0 


(15) 


Shipment  of  Units  via  Airdrop  (One  for  each  period  and 
cargo  size) 


i_|  ^ cargo ( i,k)*x(i,3,k,l) 


Y  ’ 

-  E  AB( y ) *TON ( y , k ) *u( y , 1 , 1 ) 
y=l 

P(k, 1-1 )  -  P(k, 1 )  =  0  ( 16  ) 


Fxont  Line  Supplies  (One  for  each  period) 


l  4  Y  3  1’  -My 

cargo(i,3)  *  x(i,J,3,l')  +  E  E  :  [ TRP ( y )  *  PL/ 
J=i  y=l  m=l  1=1 

Y’  3  l’-l 

^  *  E  PL*T0N(y,3)  *  u(y.m.l) 

y=l  m= 1  1=1 


2*DAPDFT]  *  u(y,m,l) 


II  MH 


I 

£  ADC (i,3)*x(i,3,3,l)-P(3,l-l)  +  P(3,l)  =  RC*PL 
i  =  l 


(  21 ) 


Using  the  assumptions  discussed  in  the  scenario 
chapter,  the  final  list  of  equations  produces  a  model  with 
168  equations  and  288  variables.  The  linear  programming 
model  was  constructed  using  Proc  IML  to  build  a  168  by  288 
matrix  consisting  of  48,384  zero  elements.  The 
approximately  2400  nonzero  elements  were  inserted  using 
matrix  notation,  and  the  matrix  was  converted  to  a  data 
set.  The  right  hand  side  values  were  constructed  as  a  two 
variable  vector  consisting  of  constraint  sign  and  right 
hand  side  value.  This  vector  was  merged  with  the  data  set 
to  form  the  objective  equation  and  the  167  constraints. 

The  linear  program  was  executed  using  the  SAS/OR 
statistical  package  (28).  SAS/OR  was  also  used  to  execute 
the  required  number  of  runs  to  build  the  dependent 
variable  outputs  required  to  construct  the  proper  seven 
variable  experimental  design  developed  by  Box  and  Behnken 
(3:461).  Once  the  required  linear  programming  runs  were 
completed,  the  SAS  regression  package  (26)  was  used  to 
develop  the  response  surface  equation.  The  response 
surface  equation  and  its  applications  are  discussed  in  the 
analysis  chapter.  SAS/Graph  was  then  used  to  construct  a 


variety  of  three-dimensional  plots  by  varying  two  of  the 
independent  variables  while  setting  the  other  variables  in 
the  model  at  a  constant  value  (27).  The  value  of  the 
objective  function  was  plotted  as  it  varied  with  the  two 
selected  independent  variables,  and  the  result  is  a  visual 
presentation  of  the  interrelationships  existing  between 
the  variables.  A  more  complete  discussion  of  the  use  of 
these  statistical  packages  and  their  applications  is 
presented  in  the  analysis  chapter. 


In  an  effort  to  evaluate  the  usefulness  of  the  model 


developed  in  the  previous  chapter,  a  scenario  centered 
around  the  Far  East  was  generated.  The  theater  of 
operation  encompasses  the  Korean  peninsula,  and  the 
military  units  were  deployed  from  the  CONUS.  More 
specifically,  Kwangju  AFB  was  the  APOD  and  Suwon  Airfield 
was  the  FOL.  The  units  were  flown  from  Travis  AFB  to  the 
Far  East.  The  route  structure  is  discussed  in  further 
detail  in  this  chapter.  The  scenario  was  greatly 
simplified,  and  the  output  of  the  model  is  meant  to 
demonstrate  the  feasibility  and  usefulness  of  the  model  as 
a  tool  in  determining  optimum  force  structuring,  not  as 
the  absolute  best  answer  to  an  actual  scenario.  Several 
simplifying  assumptions  for  the  scenario  are  listed  next. 

Assumptions 

The  assumptions  for  this  model  include: 

1)  Aircraft  will  adhere  to  standard  authorized 
ground  times  as  per  the  appropriate  MAC  55-series 
regulation . 

2)  Aircraft  are  authorized  to  concurrently  refuel 
while  uploading  or  downloading  cargo 

3)  The  enroute  times  to  transit  the  system  are  based 
upon  no  air  refueling. 


14)  The  C-141  and  C-17  are  the  only  two  aircraft 


involved  in  airdrop  operations. 

15)  The  twenty  day  period  of  deployment  can  be 
reflected  in  the  linkage  of  four  five-day  periods. 

16)  The  twenty  day  period  is  a  reasonable  time  frame 
whereby  airlift  will  have  played  its  most  critical  role. 

17)  Attrition  will  be  modeled  as  a  percentage  of 
total  missions  of  a  given  mission  and  aircraft  type. 

18)  The  Battalion  Headquarters,  F-16  squadron,  and 
Truck  Company  will  not  be  delivered  to  the  front. 

19)  All  aircraft  lost  due  to  attrition  are  assumed  to 
be  carrying  cargo. 

20)  The  C-141  fleet  is  assumed  to  possess  a  fuel 
inert  system  due  to  its  extreme  vulnerability  to  a  threat 
environment. 

21 )  OSAN  AB  is  assumed  to  be  unusable  due  to  runway 
damage . 

22)  All  bases  are  assumed  to  be  operational  24 
hours/day  * 

Input  Parameters 

The  number  of  sorties  generated  by  the  airlift  force 
depends  upon  the  number  of  aircraft  available  for  the 
operation  and  the  time  in  days  required  for  an  aircraft  to 
transit  the  system  and  return  to  home  station.  The  number 
of  aircraft  available  for  the  airlift  has  been  estimated 
.  approximate  the  normal  aircraft  by  type  available  to 

...vi  Air  Force.  It  is  assumed  that  21st  Air  Force  crews 
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and  aircraft  would  augment  22nd  Air  Force  by  flying  a 
percentage  of  the  normal  channel  missions.  This  scenario 
assumes  the  following  aircraft  will  participate  in  the 
contingency  operation  to  Korea. 

Table  I 

Aircraft  Availability 


Aircraft 

Tvoe 

Home 

Station 

Organization 

Number 

Available 

C-5A/B 

Travis  AFB 

60th  MAW 

40 

C-17 

west  coast 

? 

0 

C-141B 

McChord  AFB 

62nd  MAW 

C-141B 

Travis  AFB 

60th  MAW 

110  ( to' 

C  141B 

Norton  AFB 

63rd  MAW 

Cargo  747 

CRAF 

30 

DC-8 

CRAF 

20 

C-130 

Clark  AFB 

374th  TAW 

16 

C-130 

Yokota  AFB 

316th  TAG 

la 

C-130 

McChord  AFB 

62nd  MAW 

16 

For  the  purposes  of  the  experimental  design,  the 
number  of  C-5B’s  was  varied  from  20  to  60  available  for 
the  operation.  The  assumption  of  60  C-5’s  would  require 
21st  Air  Force  participation  in  the  airlift  from  Dover 
AFB.  The  use  of  20  or  40  C-5’s  would  assume  the  airlift 
for  this  particular  model  would  come  largely  from  Travis 
AFB,  with  the  22nd  Air  Force  C-5  crews  and  aircraft 
picking  up  the  routine  channel  missions  which  would  still 
be  required  to  be  supported.  The  number  of  C-17’s  was 
varied  from  the  existing  zero  level  to  the  addition  of  20 
or  40  aircraft.  The  C-I30’s  were  varied  from  an  average 
of  60  to  a  low  of  40  and  a  high  of  80.  The  60  and  80 
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levels  assume  additional  support  from  another  C-130  wing 
stationed  in  the  United  States,  but  placed  on  rotation  to 
the  Pacific  for  a  thirty  day  period. 

The  number  of  days  required  for  an  aircraft  to 
complete  an  entire  mission  is  dependent  upon  the  aircraft 
speed  and  the  distance  between  airfields.  The  speed  of 
the  aircraft  used  for  computations  is  not  the  true 
airspeed  (TAS)  but  rather  the  groundspeed,  which  is 
computed  for  each  mission  leg,  depending  upon  forecast 
wind  factors.  The  wind  factors  were  obtained  from  HQ 
MAC/XOSS.  Even  though  wind  factors  can  change  appreciably 
based  upon  the  time  of  year,  the  average  values  obtained 
from  MAC  will  suffice  for  the  purposes  of  the  model.  The 
intertheater  aircraft  routing  requires  an  initial 
departure  from  Travis,  an  enroute  refueling  stop  at 
Elmendorf,  and  either  a  flight  to  the  APOD  (Kwangju)  or 
direct  to  the  FOL  (Suwon) .  In  both  cases  the  crews  will 
recover  to  Yokota  Air  Base  for  refueling  and  crew  changes 
before  the  return  trip  back  through  Elmendorf  to  Travis 
AFB.  Standard  ground  times  and  an  expected  maintenance 
delay  time  were  added  to  the  figures  to  determine  an 
estimate  of  the  round  trip  time  required  for  each  aircraft 
type.  Intratheater  times  are  calculated  using  the  same 
criteria  with  the  following  results: 


Table  II 


Aircraft  Usage  Rates (days) 


Aircraft  Type 

Inter  (i) 

BBaiW 

C-5B 

2.2 

N/A 

C-17 

2.0 

0.20 

C-141 

2.1 

N/A 

C-747 

2.0 

N/A 

DC-8 

2.0 

N/A 

C-130 

N/A 

0.25 

Attrition  for  airlift  aircraft  is  a  subjective  issue 
with  little  actual  data  to  support  a  particular  viewpoint, 
yet  there  is  no  disagreement  that  the  loss  of  a 
substantial  percentage  of  our  airlift  resources  could 
prove  disastrous  to  any  major  military  operation. 

Strategic  airlift  resources  are  currently  unable  to  meet 
desired  goals  in  terms  of  metric  ton-miles  per  day  with  no 
attrition,  so  even  a  low  level  of  attrition  could  prove  to 
be  significant  to  the  war  effort.  As  pointed  out  in  a 
recent  article,  "more  than  ever  before,  aircraft  will  be 
destroyed  far  faster  than  replacements  can  be  built" 
(10:341). 
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The  value  of  determining  the  expected  attrition  of 
airlift  aircraft  is  not  so  much  to  model  the  various  usage 
rates  for  each  aircraft  type,  but  rather  to  use  aircraft 
attrition  as  an  approximation  of  the  average  loss  in 
combat  power  for  each  type  of  aircraft. 

Aircraft  attrition  has  been  approximated  by  first 
assigning  an  attrition  factor  to  each  aircraft  type  based 
upon  the  mission  to  be  flown.  The  four  mission  types  were 
assigned  a  subjective  ranking  from  zero  to  one  in  0.25 
increments  based  upon  which  mission  could  be  expected  to 
result  in  a  higher  level  of  attrition.  The  airdrop 
mission  was  assumed  to  possess  the  highest  risk  and  was 
assigned  a  value  of  one,  whereas  a  mission  to  the  APOD  was 
considered  the  safest  and  assigned  a  rating  of  0.25. 
Missions  departing  a  safe  location  and  destined  for 
another  relatively  safe  location,  such  as  the  FOL,  were 
assigned  a  value  of  0.5. 

These  rating  were  purely  subjective  and  in  reality 
enemy  forces  could  focus  on  one  particular  method  of 
delivery  and  without  near  total  air  supremacy,  cargo 
aircraft  on  any  of  these  missions  would  be  highly 
susceptible  to  enemy  fire.  Nevertheless,  these  values  do 
provide  a  basis  for  useful  comparison  for  the  purposes  of 
this  model.  Each  aircraft  type  was  also  assigned  a 
relative  ranking  based  upon  expected  survivability  in  a 
threat  environment.  The  baseline  as  far  as  the  least 
survivable  aircraft  was  established  with  the  two 


commercial  (CRAF)  aircraft,  each  receiving  a  value  of 
one.  Based  upon  such  factors  as  combat  aircrew  training, 
aircraft  size,  redundancy  of  major  systems,  and  fire 
suppressant  capability  (i.e.  foam  in  wings,  liquid 
nitrogen  systems,  etc.)  the  C-5B,  C-130,  and  C-141  were 
each  assigned  a  value  of  0.75.  An  expanded  discussion  of 
the  C-141 's  vulnerability  due  to  fuel  vapor  in  its  tanks 
will  follow  in  this  chapter.  The  C-17  was  assigned  a 
rating  of  0.5,  making  it  the  least  vulnerable  aircraft, 
due  to  proposed  design  enhancements  over  current  aircraft 
in  the  airlift  inventory.  The  two  factors  for  mission 
type  and  aircraft  type  were  combined  to  produce  an  overall 
attrition  factor  for  each  aircraft  flying  a  particular 
mission.  While  somewhat  simplistic,  this  approach 
appears  more  realistic  than  assuming  a  constant  attrition 
rate  regardless  of  aircraft  or  mission. 

As  stated  in  the  list  of  assumptions,  these  attrition 
factors  are  based  upon  the  C-141B  possessing  a  fuel  inert 
system.  The  type  of  system  is  not  at  issue,  but  while  the 
C-5B  has  a  liquid  nitrogen  system  (LN2),  and  the  C-130  has 
the  fire  suppressant  foam  in  its  wings,  the  C-141  has  no 
such  capability.  On  an  intertheater  airdrop  mission,  a 
typical  C-141B  would  probably  have  fuel  only  in  its  main 
fuel  tanks  (approximately  44,000  lbs)  by  the  time  it 
approached  the  drop  zone.  That  would  mean  the  extended 
and  aux  range  tanks  would  be  filled  with  fuel  vapor, 
waiting  for  a  source  of  ignition.  Discussions  with  Mr. 


Redman  (25)  and  Major  Maroney  (16),  both  from  ASD/XRS, 
highlighted  the  importance  of  a  fuel  inert  system  in  terms 
of  the  vulnerability  of  an  aircraft.  Their  studies 
emphasized  that  a  large  percentage  of  the  vulnerability  of 
an  aircraft  is  directly  related  to  the  inherent  risk  of 
operating  an  aircraft  in  a  threat  environment  with  fuel 
tanks  filled  with  fuel  vapor  (ullage).  The  MAC  aircraft 
most  vulnerable  to  this  ullage  risk  is  also  currently  the 
"workhorse"  of  MAC,  the  C-141B  Starlifter.  If  the  C-141B 
is  to  be  used  in  an  environment  where  attrition  is  a 
factor,  the  installation  of  a  fuel  inert  system  on  the 
C-141B  should  be  given  high  priority  by  MAC  (25).  This 
model  assumes  a  fuel  inert  system  on  the  C-141B.  The 
analysis  chapter  of  this  thesis  discusses  the  impact  in 
sorties  and  combat  power  delivered  if  the  C-141B  continued 
to  have  no  such  system. 

The  expected  survivability  of  an  aircraft  can  be 
analyzed  by  discussing  survivability  in  terms  of  its  two 
major  components,  vulnerability  and  susceptibility  (Id). 
Vulnerability  refers  to  the  characteristics  of  the 
aircraft,  such  as  redundant  flight  control  systems  (routed 
through  different  locations  of  the  fuselage),  fire 
suppressant  systems,  aircrew  combat  training,  etc. 

(19:5).  Susceptibility  refers  to  the  threat  exposure 
which  relates  more  t'  the  characteristics  of  the  mission 
(19:4).  Threats  for  a  certain  mission  type  can  be 
minimized  by  avoiding  the  threat,  or  minimizing  exposure 
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to  that  threat.  Minimizing  exposure  can  be  achieved  by 
aircraft  speed  or  perhaps  more  critically,  aircraft 
altitude  and  terrain  following  capability.  Vulnerability 
and  susceptibility  are  combined  in  an  effort  to  quantify 
survivability  of  an  aircraft.  This  simplified  aspect  of 
the  model  combines  the  susceptibility  attrition  factor 
(mission  type)  with  the  vulnerability  attrition  factor 
(aircraft  type)  to  determine  an  overall  factor  of 
attrition.  When  the  scenario  specifies  a  given  attrition 
rate,  (i.e.  five  per  cent),  each  aircraft  and  mission 
combination  is  multiplied  by  its  respective  factor  of 
attrition. 

The  number  of  aircraft  available  for  a  given  period 
is  determined  by  subtracting  the  summation  of  previous 
attrition  values  for  that  specific  aircraft  type.  The 
summation  of  attrited  aircraft  is  then  used  to  estimate 
the  loss  of  combat  power  due  to  the  loss  of  aircraft.  An 
average  value  of  combat  power  per  aircraft  load  was 
estimated  by  varying  the  number  of  aircraft  in  the  model, 
relaxing  some  of  the  constraints  to  remain  feasible,  and 
examining  the  changes  in  the  objective  function.  The 
combat  power  per  load  estimates  and  the  estimates  of  the 
number  of  loads  (aircraft)  lost  were  then  used  to 
determine  the  combat  power  lost  due  to  aircraft  attrition 
This  attrition  of  combat  power  due  to  aircraft  attrition 
was  modeled  by  including  the  variable  [ATT(i)],  the 
expected  attrition  by  aircraft  type,  in  the  objective 


function.  The  coefficient  for  each  aircraft  attrition 


variable  is  the  constant  [CPLOST(i)],  the  estimated  value 
of  combat  power  for  each  aircraft  type.  This  formulation 
allows  the  objective  function  to  be  decreased  by  the 
combat  power  lost  on  the  aircraft  which  were  attrited  by 
the  model.  The  attrition  factors  and  estimated  combat 
load  per  aircraft  estimates  are  included  in  the  following 
table : 


Table  III 

Attrition  and  Combat  Power/Load  Factors 

Estimated  Combat 


Aircraft 

tjjgajon  Attrition  Factor 

Power/Load 

C-5B 

Inter/APOD 

1.0 

0.1075 

C-17 

Inter/APOD 

0.75 

0. 131 

M 

Inter/Direct 

1.0 

" 

" 

Inter/Airdrop 

1.5 

" 

»« 

Intratheater 

1.0 

" 

C-141 

Inter/APOD 

1.0 

0.061 

’* 

Inter/Direct 

1 . 25 

" 

" 

Inter/Airdrop 

1.75 

" 

C-747 

Inter/APOD 

1 . 25 

0.143 

DC-8 

Inter/APOD 

1 . 25 

0.062 

C-130 

Intratheater 

1 . 25 

0.025 

It  may  be  noted  that  the  estimated  combat  power  per 
load  is  higher  for  the  C-17  than  it  is  for  the  C-5.  This 
factor  is  due  to  the  higher  value  CPI(l)  give'n  to  aircraft 
loads  that  arrive  directly  to  the  front.  If  the  C-5  could 
deliver  it’s  cargo  directly  to  the  objective  area,  the 
combat  power  load  factor  for  the  C-5  would  be  closer  to 
0.15.  This  factor  does  not  imply  that  all  C-17  loads 
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would  be  delivered  directly  to  the  front,  but  the  mean 
value  for  the  cargo  delivered  would  be  higher  for  the  C-17 
than  it  would  be  for  the  C-5. 

The  utilization  rate  of  the  aircraft  determines  how 
many  hours  an  aircraft  can  be  expected  to  be  airborne  in  a 
24  hour  period.  Utilization  (UTE)  rate  is  calculated  in 
peacetime  by  dividing  the  total  programmed  flying  hours  by 
the  primary  aircraft  authorized  (PAA),  then  dividing  that 
figure  by  360  to  arrive  at  a  daily  UTE  rate.  The  actual 
surge  capability  UTE  rates  are  classified,  but  they  can  be 
expected  to  exceed  current  actual  flying  hours.  Estimated 
UTE  rates  are  found  in  Table  IV.  UTE  rates  are  dependent 
upon  aircraft  type  and  type  of  mission.  In  developing  the 
model,  the  maximum  UTE  rates  varied  as  the  number  of 
aircraft  were  increased  or  decreased  in  the  model. 


Table  IV 


Aircraft  Capabilities 


J - UTE  RATES - \ 

Intertheater  Pax  Gnd  463L 


A/C  M 

30D/Direcl 

:  Intra 

Airdroo 

TAS 

Times 

Pallets 

C-5 

11/0 

0 

0 

320 

430 

(hours ) 
3.25 

36 

C-17 

14/13 

12 

10 

160 

430 

2.0 

18* 

C-141 

12.5/10 

10 

8 

155 

410 

2.25 

13** 

C-747 

10/0 

0 

0 

0 

430 

3.6 

36 

DC-8 

10/0 

0 

0 

364 

430 

2.8 

0 

C- 130 

0/0 

4 

0 

64 

270 

2.25 

6 

*  approximately  10  for  airdrop 

**  approximately  8  for  airdrop 
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The  ramp  capacity  of  an  airfield  can  play  a  critical 
role  in  determining  the  value  of  that  airfield  as  an 
operating  base.  As  determined  in  the  previous  chapter, 
ramp  saturation  occurs  when  all  available  parking  spots 
are  occupied.  The  expected  ground  time  for  each  aircraft 
helps  determine  how  many  of  that  type  aircraft  are  capable 
of  parking  on  the  airfield.  HQ  MAC/DEP  has  a  computer 
data  base  consisting  of  hundreds  of  airfields  which  have 
been  surveyed  to  determine  the  maximum  aircraft  on  the 
ground  (MOG)  based  upon  aircraft  and  airfield 
characteristics.  APOD  ramp  capacity  for  the  model  was 
varied  by  allowing  the  ramp  to  vary  in  capacity  by  25  per 
cent  from  the  actual  ramp  capacity  at  Kwangju  (i.e.  C-5 
ramp  capacity  was  varied  from  six  to  ten).  Table  V 
contains  the  ramp  capacity  by  aircraft  type. 


Table  V 


Ramp  Capacity 


Airfield 

Location 

C- 

C- 

C- 

C- 

DC- 

C- 

F- 

CAFB) 

Identifier 

££ 

17 

141B 

747 

£ 

130 

IS 

Travis 

KSUU 

42 

97 

67 

60 

60 

N/A 

N/A 

Elmendorf 

PAED 

16 

45 

39 

30 

30 

N/A 

N/A 

Yokota 

RJTY 

7 

16 

15 

8 

8 

21 

N/A 

Kwangju 

RKJJ 

8 

17 

12 

10 

10 

18 

40 

Suwon 

RKSW 

3 

15 

12 

5 

5 

17 

N/A 

80 
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Another  area  affecting  airfield  capability  is  the 
ability  to  efficiently  process  the  loading  and  unloading 
of  equipment  brought  in  by  air.  The  availability  of  heavy 
equipment,  such  as  K-loaders,  increases  the  productivity 
of  the  ALCE  personnel.  For  the  purpose  of  this  model,  the 
APOD  is  assumed  to  be  able  to  handle  1000  pallet 
equivalents  per  day,  while  the  FOL  can  handle  200  pallet 
equivalents  per  day.  Pallet  equivalents  are  used  as  a 
common  measure  of  cargo  type  since  not  all  cargo  is 
palletized.  As  the  ALCE  arrives  at  the  APOD  the 
mathematical  formulation  of  the  constraint  allows  the 
material  handling  capability  of  the  airfield  to  increase. 
The  model  varied  the  MHE  capacity  of  the  APOD  from  700  to 
1300  pallets  per  day. 

The  units  will  arrive  to  the  FOL  by  a  combination  of 
overland  delivery  from  the  APOD,  intratheater  airlift  from 
the  APOD,  intertheater  airlift  from  the  United  States,  and 
intertheater  airlift  via  airdrop  missions.  Direct 
delivery  is  provided  by  the  C-141B  and,  if  available,  the 
C-17.  Only  cargo  required  at  the  FOL  is  allowed  to  be 
shipped  via  direct  delivery.  All  unit  types  (from  1-Y) 
will  be  allowed  to  be  shipped  to  the  APOD  and  all  aircraft 
types  except  the  C-130  are  involved  in  intertheater 
airlift . 

Another  factor  varied  in  the  model  is  the  ground 
distance  from  the  APOD  to  the  FOL  expressed  in  terms  of 
days  of  travel  (in  increments  of  one  period).  While  this 
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factor  may  not  change  significantly  for  any  given 
scenario,  the  impact  of  the  distance  between  the  APOD  and 
the  FOL  is  important  and  should  be  examined.  Distance  in 
travel  periods  for  the  experimental  design  was  allowed  to 
vary  from  zero  days  between  APOD  and  FOL  to  10  days  travel 
between  the  two  fields. 

The  planned  allowable  cabin  load  (ACL)  for  an 
aircraft  depends  upon  the  type  of  cargo  to  be  carried, 
cargo  weight,  estimated  enroute  flight  time,  forecast 
weather  conditions,  and  the  availability  of  air 
refueling.  Air  refueling  would  enhance  the  payload  and 
sortie  generation  capability  of  the  C-5,  C-17,  and  C-141, 
but  for  the  purposes  of  this  model  air  refueling  is 
assumed  not  to  be  available  due  to  higher  priority 
requirements  for  the  SAC  aircraft.  The  inclusion  of  air 
refueling  for  this  scenario  would  require  recomputing  the 
enroute  time  from  Travis  AFB  to  Korea  deleting  the 
requirement  to  refuel  at  Elmendorf  AFB.  and  adding 
additional  cargo  for  a  refueling  contact  point  near 
Alaska.  The  allowable  cabin  load  for  the  aircraft  type 
was  extracted  from  data  obtained  from  the  office  of  the 
director  of  Program  Analysis  and  Evaluation,  in 
Washington,  D.C.  (9).  The  cargo  capability  in  tons  is 
found  in  the  following  table,  and  is  based  upon  a  3100 
nautical  mile  critical  leg. 
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Table  VI 


Cargo 

Capability 

(tons ) 

Aircraft  Tvoe 

Outsize 

Oversize 

Bulk 

Personnel 

(people ) 

C-5B 

53 

59 

83 

320 

C-17 

47 

32 

37 

160 

C-141 

00 

25 

30 

155 

C-747 

00 

71 

97 

00 

DC-8 

00 

00 

41 

364 

C-130 

00 

00 

13.8 

64 

Units  required  to  be  airlifted  are  dependent  upon  the 
scenario.  This  model  assumed  the  same  units  as  described 
in  Cooke’s  thesis  (7:83).  The  following  data  was  extracted 
from  FM  101-10-1  and  Cooke’s  thesis: 


- 1 
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THE  USE  OF  HATHEHATICAL  PROGRAHHING  AND  RESPOOSE 
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Table  VII 


Units  To  Deploy  (tons) 


Tyne  Number 

Outsize 

Airborne 

Battalion 

9 

0 

Airborne 
Battalion  HQ 

3 

73 

Air  Assault 
Battalion 

9 

152 

Mechanized 

Battalion 

9 

2055 

155mm  Art. 
Battalion 

3 

139 

F-16  Sqdn  * 

3 

0 

Medium  Truck 
Company  ** 

3 

1098 

ALCE 

2 

176 

*  18  aircraft 
**  180,000  ton-miles/day 

Daily 

Supply 


Oversize 

Bulk 

Pax 

Read 

1400 

268.7 

1400 

149.2 

0 

0 

305 

0 

990.3 

574.7 

1960 

446 

3885 

151.7 

2052 

400.2 

1243 

91.5 

710 

212.85 

249.6 

155.4 

472 

102.5 

405 

270.0 

250 

55.7 

2213 

244.0 

400 

73.8 

In  order  to  quantify  the  amount  of  combat  power 
delivered  to  the  front,  the  following  inputs  must  be 
examined:  the  capability  of  the  aircraft  to  deliver  the 
cargo,  the  quantity  of  men  and  equipment  per  unit  to  be 
delivered,  the  values  attached  to  the  three  MOE’s  of 
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combat  power  and  the  scaling  developed  for  delivering 
combat  power  over  time.  The  first  two  have  been 
discussed,  and  the  next  two  were  developed  from  previous 
thesis  efforts  and  the  Army  Command  and  General  Staff 
College  (7:83).  The  data  is  presented  in  the  following 
tables : 

Table  VIII 

Measures  of  Effectiveness 


Type 

Fire 

P0£££ 

Anti-Tank 

Strength 

Defensive 

Frontage 

Airborne  Battalion 

4 

19.5 

4 

Airborne  Batt .  HQ 

0 

0 

0 

Air  Assault  Batt. 

6 

28.5 

4 

Mechanized  Batt. 

8 

40 

6 

155  Art.  Battalion 

3 

3 

0 

F-16  Sqdn 

8 

36 

0 

Medium  Truck  Company  0 

0 

0 

ALCE 

0 

0 

0 

Table  IX 


Combat  Value  Over  Time 


Period  Combat  Value  Factor 


1 

2 

3 

4 


2.5 

1.8 

1.3 

1.1 


Units  destined  for  the  front  that  are  airdrop  capable 
are  included  in  the  scenario,  along  with  the  only  two 
aircraft  capable  of  intertheater  airdrop  missions,  the 
C-141B  and  the  proposed  C-17.  The  scenario  allows  the 
C-141  to  airdrop  supplies  and  personnel,  while  the  C-17  is 
assumed  to  possess  the  capability  to  airdrop  all  four 
types  of  cargo.  While  not  yet  operational,  the  C-17  is 
being  developed  with  the  capability  to  airdrop  outsized 
cargo,  and  that  capability  is  included  in  the  development 
of  this  scenario  (2).  The  quantity  of  cargo  delivered  by 
airdrop  operations  is  constrained  by  the  total  number  of 
airdrop  capable  aircraft  available,  the  number  of  pallets 
per  aircraft  in  the  airdrop  configuration,  and  the  rigging 
capability  of  the  army  riggers.  This  model  assumes  an 
army  cargo  rigging  capability  of  400  pallets  per  day,  and 
it  also  assumes  that  all  C-141B  and  C-17  are  airdrop 
capable  aircraft.  Table  IV  lists  the  pertinent  airdrop 


data. 


Having  determined  the  input  parameters  through  a 
combination  of  historical  data  and  a  subjective  estimation 
of  the  range  of  values  for  the  scenario  generated  for  this 
model,  the  following  chapter  will  analyse  the  linear 
program  and  response  surface  outputs. 
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V.  Analysis 


Introduction 

While  linear  programming  remains  a  powerful  and  viable 
analytical  tool,  it  does  have  limitations.  The  standard 
sensitivity  analysis  of  ranging  the  right  hand  sides  allows 
for  only  single  parameter  analysis.  The  restriction  of 
varying  only  one  parameter  at  a  time  would  require  an 
inordinate  number  of  model  runs  to  fully  examine  the 
sensitivity  of  the  model  to  changing  inputs.  An 
alternative  methodology,  developed  by  Smith  and  Mellichamp 
(29),  allows  the  use  of  response  surface  methodology  and 
multivariate  regression  analysis  to  perform  a 
multidimensional  impact  analysis  of  the  model. 

This  methodology  uses  regression  analysis  as  a  curve 
fitting  tool  to  estimate  the  coefficients  of  the 
mathematical  equation  used  to  describe  the  response 
surface.  The  generated  response  surface  provides  a  means 
of  depicting  the  different  relationships  between  the 
factors  of  the  model. 

Experimental , Design  Development 

The  experimental  design  developed  by  Box  and  Behnken 
(3:461)  was  selected  due  to  its  rotatability  and  nearly 


orthogonal  design.  Even  though  the  constant  term  and 


the  quadratic  estimates  0^  are  correlated,  the 


correlation  is  negligible  (see  Table  X),  and  the  overall 
design  still  possesses  a  high  degree  of  orthogonality.  The 
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three  level  rotatable  design  included  seven  variables 
requiring  only  56  points  plus  an  added  point  at  the 
origin.  This  design  provided  efficient  estimates  of  the 
P’s  for  the  second  order  polynomial. 

The  design  was  generated  by  selecting  a  range  for  the 
seven  variables  to  be  investigated,  coding  the  variables, 
and  running  the  linear  program  57  times  to  generate  the 
required  objective  function  values.  The  notation  for  this 
particular  design  requires  the  lower  limit  of  the 
individual  variable  to  be  coded  as  a  "-1" ,  the  center  point 
value  of  the  variable  to  be  coded  as  a  "0",  and  the  upper 
limit  value  to  be  coded  with  a  "+1" .  The  coded  design  is 
depicted  in  Appendix  A.  The  seven  variables  selected  to  be 
investigated  were:  1)  the  number  of  C-5B  aircraft 
available,  2)  the  number  of  017  aircraft  available,  3) 

APOD  MHE  capability,  4)  aircraft  attrition,  5)  distance 
from  the  APOD  to  the  FOL,  6)  APOD  ramp  space,  and  7)  the 
number  of  0130’s  available  for  the  operation.  The 
noncoded  variable  values  are  listed  in  Appendix  B.  The 
interaction  and  quadratic  terms  were  generated  using  SAS. 
The  SAS  program  is  listed  in  Appendix  C.  The  interaction 
terms,  such  as  xl2,  refer  to  the  interaction  between  two 
particular  main  effects.  For  this  example,  the  interaction 
is  between  05’s(xl)  and  017’s(x2).  In  three  dimensions, 
the  plane  generated  by  the  xl  and  x2  variables  is  twisted 
by  the  interaction  term  xl2.  The  magnitude  of  the  twist  is 
determined  by  the  coefficient  for  x!2.(7:99) 


Quadratic  terms,  such  as  xll,  refer  to  the  interaction  of  a 
main  effect  term  with  itself.  The  coefficient  of  a 
quadratic  response  expresses  the  degree  of  curvature  of  the 
response  surface  to  that  particular  axis.  For  the  full 
model,  there  are  seven  main  effect  terms,  21  interaction 
terms,  seven  quadratic  terms,  and  the  intercept  term,  for  a 
total  of  36  terms. 

Response  Surface  Development 

Utilizing  a  second  order  model  with  the  appropriate 
interaction  and  quadratic  terms ,  a  response  surface  was 
generated.  The  linear  programming  model,  consisting  of  168 
constraints  and  288  variables,  was  run  57  times  using  the 
points  of  the  design  depicted  in  Appendix  B.  The  seven 
variables  were  abbreviated  using  the  following  format: 


Variables 

Range 

xl 

- 

Number  of  C-5B  aircraft 

20 

-  60  aircraft 

x2 

s 

Number  of  C-17  aircraft 

0 

-  40  aircraft 

x3 

- 

APOD  MHE  capability 

700 

-  1300  pallets/day 

x4 

= 

Aircraft  Attrition 

0 

-  10  per  cent 

x5 

= 

Distance  APOD  to  FOL 

0 

2  periods 

x6 

= 

APOD  ramp  size  capability 

0.75 

-  1.25  std  ramp 

x7 

- 

Number  of  C-130  aircraft 

40 

-  80  aircraft 

The  interaction  and  quadratic  terms  are  as  specified 


in  the  SAS  program  in  Appendix  C. 


Using  Proc  Stepwise  with  the  "maxr"  option,  the  SAS 
program  brought  each  of  the  variables  into  the  predicted 
equation.  The  program  determined  the  best  model  based  upon 
the  number  of  variables  allowed  in  the  program.  Since 
there  is  no  variance  in  the  output  of  a  deterministic 
model,  all  the  terms  were  initially  included  in  the  model. 
Proc  Stepwise  then  determined  the  best  one  variable  model, 
the  best  two  variable  model,  and  so  on  until  all  variables 
were  included.  A  regression  was  run  with  all  the  variables 
in  the  model  with  the  results  listed  in  Figure  1 .  The 
prediction  capability  of  the  full  model  is  shown  in  Table 
XI,  which  lists  the  linear  programming  output  "y",  the 
prediction  value  of  the  experimental  design,  and  the 
percentage  error  for  each  of  the  57  runs .  From  a  combat 
value  range  of  113.05  to  297.63,  the  full  model  had  a  mean 
percentage  error  of  2.247  percent. 

For  comparison,  a  model  with  less  than  the  total 

2 

variables  was  analyzed.  Using  R  ,  which  measures  the 
percentage  of  the  variation  of  the  dependent  variable  "y" 
which  can  be  explained  by  the  independent  variables 
currently  in  the  model,  and  the  error  sum  of  squares  (SSE), 
a  model  with  12  variables  was  selected  for  comparison. 

Proc  Stepwise,  with  the  "maxr”  option,  selected  the  best  12 
variable  model.  Proc  Rsquare  was  also  used  to  confirm  that 
the  correct  variables  were  selected  for  the  model.  A 
regression  of  the  12  variable  model  produced  the  results 
shown  in  Figure  2.  Table  XII  lists  the  dependent  variable 


SAS  1 

i 6 1  S2  TUESDAY,  NOVEMBER  11,  1986 


DEP  VARIABLEi 

Y 

COMBAT  FIRE  POWER 

ANALYSIS  OF 

VARIANCE 

SUM  OF 

MEAN 

SOURCE  DF 

SQUARES  SQUARE  F  VALUE 

PRQB>F 

MODEL  35 

111141.3  3175.465  35.492 

0.0001 

ERROR  21 

1878.869  89. 

46996 

C  TOTAL  56 

113020.1 

ROOT  MSE 

9.458856  R-SQUARE  0.9834 

DEP  MEAN 

196.6556  ADJ 

R-SQ  0.9557 

C.V. 

4.B09859 

PARAMETER 

ESTIMATES 

PARAMETER 

STANDARD 

T  FOR  HOi 

VARIABLE 

DF 

ESTIMATE 

ERROR  1 

9ARAMETER«0 

INTERCEP 

1 

194.185 

9.458856 

20.529 

XI 

1 

0.7795417 

1.930781 

0.404 

X2 

1 

22.25063 

1.930781 

11.524 

X3 

1 

12.78904 

1.930781 

6.624 

X4 

1 

-38.0216 

1.930781 

-19.692 

X5 

1 

-43.9425 

1.930781 

-22.759 

X6 

1 

1.234375 

1.930781 

0.639 

n 

1 

-0.00104167 

1.930781 

-0.001 

XU 

1 

-3.06 

3.861562 

-0.792 

X12 

1 

-0.69625 

3.344211 

-0.208 

X 13 

1 

-2.447 

3.344211 

-0.732 

X14 

1 

1.1555 

3.344211 

0.346 

X15 

1 

-5.95075 

3.344211 

-1.779 

X16 

1 

-0.150875 

3.344211 

-0.045 

X 17 

1 

-0.003125 

3.344211 

-0.001 

X22 

1 

-14.0293 

3.861562 

-3.633 

X23 

1 

2.355625 

3.344211 

0.704 

X24 

1 

-6.715 

3.344211 

-2.008 

X  25 

1 

12.361 

3.344211 

3.696 

X26 

1 

-0.526875 

3.344211 

-0. 158 

X27 

1 

0 

3.344211 

0.000 

X33 

1 

-3.20913 

3.861562 

-0.831 

X34 

1 

-3.71075 

3.344211 

-1.110 

X35 

1 

-8.7935 

3.344211 

-2.629 

X36 

l 

1.569875 

3.344211 

0.469 

X37 

1 

0 

3.344211 

0.000 

X44 

l 

3.7145 

3.861562 

0.962 

X45 

1 

0.036625 

3.344211 

0.011 

X46 

l 

0. 101875 

3.344211 

0.030 

X47 

1 

0 

3.344211 

0.000 

X55 

1 

22.21737 

3.861562 

5.753 

X56 

1 

-0.071625 

3.344211 

-0.021 

X57 

1 

0 

3.344211 

0.000 

X66 

1 

-0.44325 

3.861562 

-0.115 

X67 

1 

-0.003125 

3.344211 

-0.001 

X77 

1 

0.677375 

3.861562 

0.  175 

Fig  1.  Rtgrtttion  Output  -  Full  Modil 


SAS  8 

1 6 i 52  TUESDAY,  NOVEMBER  il,  1986 

DEP  VARIABLE]  Y  COMBAT  FIRE  POWER 

ANALYSIS  OF  VARIANCE 

SUM  OF  MEAN 

SOURCE  OF  SfiUARES  SQUARE  F  VALUE  PRQB>F 

MODEL  12  110762.6  9230.217  179.900  0.0001 

ERROR  44  2257.323  51.30739 

C  TOTAL  56  113020.1 

ROOT  MSE  7.162918  R-SQUARE  0.9800 

DEP  MEAN  196.6556  ADJ  R-SQ  0.9746 

C.V.  3.642367 

PARAMETER  ESTIMATES 


VARIABLE 

DF 

PARAMETER 

ESTIMATE 

STANDARD 

ERROR 

T  FOR  HOi 
PARAMETER-0 

INTERCEP 

l 

189.9867 

1.928193 

98.531 

X2 

1 

22.25063 

1.462125 

13.218 

X3 

1 

12.78904 

1.462125 

8.747 

X4 

1 

-38.0216 

1.462125 

-26.004 

X5 

1 

-43.9425 

1.462123 

-30.054 

X13 

1 

-3.93073 

2.332474 

-2.350 

X22 

l 

-12.7173 

1.97381 

-6.436 

X24 

1 

-6.713 

2.532474 

-2.652 

X25 

1 

12.361 

2.532474 

4.881 

X34 

1 

-3.71075 

2.532474 

-1.463 

X35 

1 

-8,7935 

2.532474 

-3.472 

X44 

1 

5.026457 

1.97581 

2.544 

XS5 

1 

23.32933 

1.97581 

11.909 

Fig  2.  Rtgrmion  Output  -  Sfltcttd  Modtl 


Tab  1 


X  I 


SAS  S 

1 6 1  52  TUESDAY |  NOVEMBER  119  1986 


OBS 

ERR0R1 

Y 

FRED 

1 

0.036797 

144.126 

139.011 

2 

0.038448 

141.729 

136.482 

3 

-0.044587 

216.846 

226.966 

4 

-0.040437 

215.086 

224. 150 

5 

0.042202 

223.841 

214.777 

6 

0.047387 

222.775 

212.655 

7 

-0.017325 

297.631 

302.878 

8 

-0.017024 

293.335 

300.470 

9 

0.005771 

194.330 

193.215 

10 

0.005825 

194.355 

193.229 

11 

0.003192 

191.664 

191.054 

12 

0.003181 

191.664 

191.056 

13 

-0.005894 

188.074 

189. 189 

14 

-0.003949 

188.074 

189.200 

13 

-0.003166 

191.356 

191.964 

16 

-0.003177 

191.336 

191.966 

17 

-0.044392 

183.119 

193.719 

18 

-0.044403 

185.119 

193.721 

19 

0.044535 

268.329 

256.884 

20 

0.044363 

268.329 

256.882 

21 

0.036278 

243.704 

237. 102 

22 

0.036269 

245.704 

237.104 

23 

-0.091949 

113.050 

124.497 

24 

-0.091934 

113.030 

124.493 

23 

0.032337 

164.738 

159.563 

26 

-0.001673 

246.312 

246.725 

27 

-0.008920 

188.490 

190.187 

28 

0.037932 

134.811 

129.884 

29 

0.003316 

123.035 

124.622 

30 

-0.027407 

184.351 

189.546 

31 

-0.019796 

243.943 

248.870 

32 

0.010800 

138.782 

157.083 

33 

-0.000195 

166.391 

166.423 

34 

-0.000207 

166.391 

166.423 

33 

0.006686 

251.561 

249.890 

36 

0.006694 

251.561 

249.888 

37 

0.000139 

216.923 

216.889 

38 

0.000149 

216.923 

216.891 

39 

-0.011283 

146.596 

148.269 

40 

-0.011269 

146.596 

148.267 

41 

-0.048024 

154.761 

162.568 

42 

-0.004392 

278.712 

279.942 

43 

-0.040466 

231.891 

241.670 

44 

0.046063 

166.817 

159.471 

43 

0.007455 

166. 149 

164.919 

46 

0.034960 

231.123 

223.316 

47 

-0.027069 

264.029 

271.375 

48 

0.055001 

187.584 

177.805 

49 

-0.003682 

215.380 

216. 176 

30 

-0.003492 

210.459 

211.621 

31 

-0.013778 

181.930 

184.472 

32 

-0.013010 

180.004 

182.747 

S3 

0.008071 

145. 173 

144.011 

95 


Tablt  XI (continutd) 


SAS  6 

1 6 i 52  TUESDAY,  NOVEMBER  11,  1986 


OBS 

ERR0R1 

Y 

PRED 

34 

0.0055424 

144.424 

143.628 

55 

0.0169993 

164.098 

161.355 

36 

0.0151277 

170.539 

168.017 

37 

-0.0000000 

194.183 

194.185 

VARIABLE 


Y 

FRED 

ERR1 


i 

H 


i 


LABEL 

COMBAT  FIRE  POWER 
PREDICTED  VALUE 


SAS  7 

16i 32  TUESDAY,  NOVEMBER  11,  19S6 


MEAN 


196.6335789 
196. 65557B9 
0.0224733 


MINIMUM 

VALUE 

113.030000 

124.493375 

0.000000 


MAXIMUM 

VALUE 

297.631000 

302.878438 

0.091949 


"y"  from  the  linear  program,  the  prediction  value  based 
upon  the  selected  model  and  the  percentage  error  for  the  57 
runs.  Using  only  12  of  the  35  variables,  the  model 
produced  a  mean  percentage  error  rate  of  only  2.687 
percent,  only  0.4  percentage  error  increase  over  the  full 
35  variable  model  (Figure  1).  By  comparison,  using  linear 
programming  to  examine  the  effects  of  the  seven  variables 
by  varying  the  aircraft  in  increments  of  five,  the  APOD  MHE 
in  increments  of  100  pallets/day,  attrition  by  .5  percent, 
ramp  space  by  0.25  percent,  and  APOD  to  FOL  distance  by  1 
day  increments  would  require  1,843,200  separate  runs  of  the 
linear  programming  model. 

Verification 

The  12  variable  model  includes  the  main  effects  of: 

1)  number  of  C-17’s,  2)  APOD  MHE  capability,  3)  aircraft 
attrition,  and  4)  the  number  of  days  travel  between  the 
APOD  and  FOL.  As  expected,  the  P  coefficients  for  the 
number  of  C-17’s  and  the  amount  of  MHE  capability  have 
positive  signs  indicating  their  positive  effect  on  the 
objective  function.  Attrition  and  the  distance  from  the 
APOD  to  the  FOL  have  negative  signs  indicating  their 
negative  effect  upon  the  delivery  of  combat  power.  As  an 
additional  check  of  the  proposed  model ,  the  percentage 
error  was  plotted  against  the  predicted  value  with  the 
results  listed  in  Figure  4  for  the  selected  model  and 
Figure  3  for  the  full  model.  The  plots  are  quite  good  with 
only  one  or  two  points  approaching  the  eight  percent  error 
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XII 


T  *  to  1  • 

SAS  12 

1 6 i S3  TUESDAY,  NOVEMBER  11,  1986 

E 


R 


H 

0  R 

X 

X 

X 

X 

X 

X 

X 

X 

P 

B  0 

X 

X 

X 

X 

1 

2 

2 

2 

3 

3 

4 

3 

R 

S  R 

2 

3 

4 

5 

3 

2 

4 

5 

4 

5 

4 

5 

Y 

E 

1  0.0SS262 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

1 

1 

144.126 

136.578 

2  0.037711 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

l 

1 

141.729 

136.578 

i 

3  -0.033936 

0 

0 

1 

-1 

0 

0 

0 

0 

0 

0 

1 

t 

216.846 

224.463 

i 

4  -0.041777 

0 

0 

1 

-1 

0 

0 

0 

0 

0 

0 

1 

1 

215.086 

224.463 

i 

3  0.052767 

0 

0 

-1 

1 

0 

0 

0 

0 

0 

0 

1 

1 

223.841 

212.622 

6  0.047733 

0 

0 

-1 

1 

0 

0 

0 

0 

0 

0 

1 

1 

222.775 

212.622 

i 

7  -0.009369 

0 

0 

-1 

-1 

0 

0 

0 

0 

0 

0 

l 

1 

297.631 

300.507 

5 

8  -0.017143 

0 

0 

-1 

-1 

0 

0 

0 

0 

0 

0 

1 

1 

295.355 

300.307 

ft 

9  0.022861 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

194.330 

189.987 

; 

10  0.022992 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

194.335 

189.987 

:s 

11  0.008828 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

191.664 

189.987 

: 

12  0.008828 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

191.664 

1B9.987 

E 

13  -0.010068 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

188.074 

189.987 

14  -0.010068 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

188.074 

189.987 

R 

13  0.007207 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

191.356 

189.987 

ft 

16  0.007207 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

191.356 

189.987 

? 

17  -0.033139 

1 

0 

0 

1 

0 

1 

0 

1 

0 

0 

0 

1 

185.119 

191.468 

% 

18  -0.033139 

1 

0 

0 

1 

0 

1 

0 

1 

0 

0 

0 

1 

185.119 

191.468 

■t 

19  0.053796 

1 

0 

0 

-1 

0 

1 

0 

-1 

0 

0 

0 

1 

268.329 

254.631 

i. 

20  0.033796 

1 

0 

0 

-1 

0 

1 

0 

-1 

0 

0 

0 

1 

268.329 

254.631 

•l 

21  0.046210 

-1 

0 

0 

-1 

0 

1 

0 

1 

0 

0 

0 

1 

243.704 

234.852 

22  0.046210 

•  1 

0 

0 

-1 

0 

1 

0 

1 

0 

0 

0 

1 

245.704 

234.852 

23  -0.073213 

-1 

0 

0 

1 

0 

t 

0 

-1 

0 

0 

0 

1 

113.030 

122.245 

& 

24  -0.073213 

•  i 

0 

0 

1 

0 

1 

0 

-1 

0 

0 

0 

1 

113.050 

122.243 

25  0.030962 

1 

0 

1 

0 

0 

1 

1 

0 

0 

0 

1 

0 

164.738 

139.810 

; 

26  -0.011919 

1 

0 

•*  1 

0 

0 

1 

-1 

0 

0 

0 

1 

0 

246.312 

249.283 

V 

27  -0.014936 

-l 

0 

-1 

0 

0 

1 

1 

0 

0 

0 

l 

0 

188.490 

191.332 

V 

28  0.047168 

-1 

0 

1 

0 

0 

1 

-1 

0 

0 

0 

1 

0 

134.811 

128.739 

D 

29  -0.028769 

-1 

0 

1 

0 

0 

1 

-l 

0 

0 

0 

1 

0 

125.035 

128.739 

30  -0.036387 

-1 

0 

0 

0 

1 

1 

0 

0 

0 

1 

0 

184.351 

191.332 

31  -0.021422 

l 

0 

-l 

0 

0 

1 

-1 

0 

0 

0 

1 

0 

243.943 

249.283 

V 

32  -0.006432 

1 

0 

1 

0 

0 

1 

1 

0 

0 

0 

1 

0 

158.782 

159.810 

• 

33  0.001934 

0 

1 

1 

0 

0 

0 

0 

0 

l 

0 

1 

0 

166.391 

166.070 

34  0.001934 

0 

1 

1 

0 

0 

0 

0 

0 

1 

0 

1 

0 

166.391 

166.070 

35  0.008121 

0 

1 

-1 

0 

0 

0 

0 

0 

-1 

0 

1 

0 

251.561 

249.535 

36  0.008121 

0 

1 

-1 

0 

0 

0 

0 

0 

-1 

0 

1 

0 

251.561 

249.533 

! 

37  0.001792 

0 

-1 

“1 

0 

0 

0 

0 

0 

1 

0 

1 

0 

216.923 

216.535 

38  0.001792 

0 

- 1 

-1 

0 

0 

0 

0 

0 

1 

0 

1 

0 

216.923 

216.535 

39  -0.008906 

0 

-1 

1 

0 

0 

0 

0 

0 

-1 

0 

1 

0 

146.596 

147.913 

.* 

40  -0.008906 

0 

-1 

1 

0 

0 

0 

0 

0 

-1 

0 

1 

0 

146.596 

147.913 

41  -0.076706 

0 

1 

0 

1 

1 

0 

0 

0 

0 

l 

0 

1 

154.761 

167.618 

f 

t 

42  -0.022033 

0 

1 

0 

-1 

- 1 

0 

0 

0 

0 

-1 

0 

1 

278.712 

284.992 

\k 

43  -0.041086 

0 

-l 

0 

-l 

-l 

0 

0 

0 

0 

1 

0 

1 

231.891 

241.827 

44  0.043040 

0 

-1 

0 

1 

1 

0 

0 

0 

0 

-1 

0 

1 

166.817 

159.627 

43  -0.031364 

0 

- 1 

0 

l 

-l 

0 

0 

0 

0 

-1 

0 

l 

166.149 

171.329 

46  0.003209 

0 

- 1 

0 

-1 

1 

0 

0 

0 

0 

1 

0 

1 

231.123 

229.925 

47  -0.033181 

0 

1 

0 

-l 

1 

0 

0 

0 

0 

-1 

0 

l 

264.029 

273.090 

48  0.044920 

0 

1 

0 

1 

- 1 

0 

0 

0 

0 

1 

0 

1 

187.584 

179.520 

! 

49  0.014464 

1 

1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

215.380 

212.309 
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E 

R 

R 


0 

R 

X 

X 

X 

X 

X 

X 

X 

X 

P 

B 

0 

X 

X 

X 

X 

1 

2 

2 

2 

3 

3 

4 

3 

R 

S 

R 

2 

3 

4 

3 

5 

2 

4 

3 

4 

5 

4 

3 

Y 

E 

SO 

-0.008714 

1 

1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

210.439 

212.309 

SI 

-0.025711 

1 

-1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

181.930 

186.731 

32 

-0.03602S 

1 

-l 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

180.004 

186.731 

S3 

0.020693 

-1 

-1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

143.173 

142.230 

34 

0.013427 

-1 

-1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

144.424 

142.230 

SS 

-0.022108 

-l 

1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

164.098 

167.808 

36 

0.016395 

-1 

1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

170.559 

167.808 

57 

0.022098 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

194.185 

189. 9B7 

13 

1986 
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VARIABLE 

LABEL 

MEAN 

MINIMUM 

MAXIMUM 

VALUE 

VALUE 

X2 

C-17  AIRCRAFT 

0.000000 

1.000000 

X3 

APOD  MHE 

0.000000 

-1.000000 

1.000000 

X4 

AIRCRAFT  ATTRITION 

1.000000 

X3 

DI8T  TO  FOL 

0.000000 

-l. 000000 

1.000000 

X1S 

1.000000 

X22 

0.421033 

rK7v  i  a  i  i  3 

1.000000 

X24 

11^1*14*1*1*  *  '!*■ 

SxJEi  •  •  In 

1.000000 

X25 

■Mil] 

-1.000000 

1.000000 

X34 

■B  ■■ 

1.000000 

X35 

0.000000 

-1.000000 

1.000000 

X44 

0.421033 

M  H  II II 1  B 

1.000000 

X55 

0.421053 

0.000000 

1.000000 

Y 

COMBAT  FIRE  PONER 

196.633579 

113.050000 

297.631000 

PRE 

PREDICTED  VALUE 

196.653379 

122.244694 

300.506611 

ERR 

0.026871 

0.001792 

0.076706 

* 

i 

* 

i 

I 

i 

( 

i 

i 

i 

I 


ERROR  I 
0.08  ♦ 


PLOT  OF  ERRQRlaPRED 


SYMBOL  USED  IS  E 


0.06  ♦ 


0.04  ♦ 


0.02  ♦ 


0.00  ♦ 


-0.02  ♦ 


-0.04  ♦ 


-0.06  ♦ 


-0.08  ♦ 


-0.10  * 


E 

E  E 


E 

E 


E  E 


E  E 


E 

E 


E 

E 


EE 
— E* 
EE 
E 


E  E 


EE 


— + — - — ♦ - ♦ - ♦ - + - ♦ - 

120  140  160  180  200  220  240 
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Fig  3.  Pvrcfntagt  Error  vi  Prtdicttd  Valuv-Full  Nodal 
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rate  for  the  selected  model.  The  combination  cf  40  C-5’s, 

0  C-17’s,  1000  pallets/day  for  MHE,  five  percent  attrition, 
and  10  days  travel  between  APOD  and  FOL,  produces  a 
response  surface  output  approximately  eight  percent  higher 
than  the  output  of  the  linear  programming  model.  Overall, 
the  percentage  error  rate  is  approximately  2.68  percent, 
which  means  the  response  surface  is  quite  accurate  for 
prediction  purposes .  Approximately  98  percent  of  the 
variability  in  the  objective  function  value  can  be 
explained  by  the  12  terms  currently  in  the  selected  model. 
In  addition,  to  the  predicted  outputs  generated  by  the  57 
runs  of  the  LP  model,  ten  mid-range  and  extreme  points  were 
generated  for  verification  purposes.  Even  with  these 
points,  the  mean  percentage  error  rate  was  only  3.68 
percent.  The  ten  runs,  model  output,  and  predicted  output, 
are  shown  in  Table  XIII. 

Table  XIII 

Mid-range  and  Extreme  Point  Data  * 


c- 

C- 

APOD 

Attri¬ 

APOD 

C- 

£ls 

IZla 

MHE 

tion 

Dist 

Bamp 

130’ 

s  Y. 

prqd 

30 

30 

1100 

0.05 

1 

1.00 

60 

201.14 

202.19 

50 

30 

1000 

0.05 

2 

1 . 25 

40 

183.42 

195.60 

60 

5 

1300 

0.10 

0 

1.25 

80 

227 . 56 

220.90 

25 

40 

800 

0.00 

0 

0.75 

40 

285.83 

275.63 

35 

15 

900 

0.05 

1 

1.10 

80 

185.34 

179.37 

35 

35 

1250 

0.00 

1 

0.75 

40 

252.26 

260.84 

50 

20 

700 

0.10 

1 

1.00 

50 

147.37 

147.91 

45 

35 

1000 

0.10 

2 

0.75 

60 

148.59 

‘  162.25 

20 

25 

1200 

0.00 

1 

1.00 

80 

244.41 

250.48 

30 

10 

1200 

0.05 

1 

1.25 

40 

186.16 

175.68 

*  Average  percentage  error  =  3.68  percent 


With  an  average  error  rate  of  only  3.68  percent  for 
the  mid-range  and  extreme  point  values  and  a  2.68  percent 


error  rate  for  the  57  model  runs,  this  response  surface 
captures  over  98  percent  of  the  variability  in  the  model 
and  verifies  that  this  quadratic  response  surface  is  an 
accurate  indicator  of  the  linear  programming  model 
responses.  With  over  a  96  percent  level  of  accuracy  over  a 
model  with  more  than  1,843,200  possible  points,  the 
advantages  of  using  the  response  surface  equation  to 
approximate  the  deterministic  model  outputs  are  obvious. 

While  the  C-5  variables  were  included  in  the 
regression  model,  only  the  C-5  interaction  with  APOD  to  FOL 
distance  was  selected  as  a  significant  C-5  variable  in 
determining  the  value  of  the  objective  function.  This  is 
not  meant  to  imply  that  the  C-5  is  not  a  significant 
asset.  Within  this  model,  the  number  of  C-5  aircraft 
available  ranges  between  20  and  60,  and  this  availability 
of  C-5’s  is  sufficient  to  preclude  it  from  becoming  a  major 
factor  is  estimating  the  value  of  the  objective  function. 

Validation  of  the  linear  program  and  the  response 
surface  is  a  more  difficult  issue.  There  has  been  no  real 
world  event  similar  to  the  scenario  generated  for  this 
thesis.  The  C-17  is  not  yet  operational  so  its  impact  upon 
a  given  scenario  can  only  be  hypothesized.  Certain  aspects 
of  the  scenario  parallel  results  of  previous 
contingencies.  The  distance  between  the  APOD  and  FOL  and 


any  aircraft  attrition  adversely  affects  the  timely 
delivery  of  combat  power.  In  over  90  percent  of  the  57 
model  runs,  MHE  capacity  at  the  APOD  was  a  binding 
constraint,  which  reflects  the  importance  of  adequate  MHE. 
This  reality  was  highlighted  in  the  1981  Air  Force  MHE 
Conference,  which  stated  that  MAC  was  significantly  short 
of  its  required  MHE  capability  required  to  support  a 
wartime  scenario .( 17 : 62 )  Throughout  most  of  the  various 
ranges  of  the  model,  the  C-141  sorties  and  usage  rates  were 
shown  to  be  binding  constraints.  Due  to  its  flexibility, 
range,  and  cargo  capacity,  the  C-141  has  been  designated  as 
the  "workhorse"  of  MAC,  and  the  model  validates  that  point 
of  view.  The  C-141  was  not  selected  as  a  critical 
variable  required  to  explain  the  variance  in  the  model 
because  of  the  range  of  aircraft  expected  to  be  available 
for  future  operations .  Regardless  of  whether  the  aircraft 
are  operated  by  active  duty,  reserve,  or  Air  Guard  crews, 
the  C-141  will  remain  in  the  airlift  inventory  for  many 
years  to  come.  This  does  not  mean  the  C-141  is  not 
important,  only  that  it  is  expected  to  be  available  in 
sufficient  numbers  that  it’s  availability  has  been 
determined  not  to  be  an  issue  for  this  particular  model. 

A  general  comparison  can  be  made  with  Cooke’s  thesis, 
but  so  many  variables  have  changed  that  a  direct  comparison 
would  be  inaccurate.  Factors  such  as  aircraft 
availability,  constraint  formulations,  UTE  rates,  the 
reformulation  of  the  excess  supply  variables,  the  inclusion 


of  aircraft  attrition,  the  inclusion  of  combat  power 
degradation  due  to  aircraft  attrition,  different  allowable 
cabin  loads  (ACL’s),  and  a  heavier  emphasis  on  goal 
programming  by  Cooke  make  a  direct  comparison  unrealistic. 
In  general,  both  efforts  provide  similar  insight  into  the 
importance  of  MHE  and  the  probable  impact  of  the  C-17  on 
the  timely  delivery  of  combat  power  to  the  theater 
commanders.  Additionally,  this  thesis  highlights  the 
critical  impact  of  distance  between  the  APOD  and  the  FOL, 
and  the  affect  of  aircraft  attrition  on  the  timely  delivery 
of  combat  power.  Several  other  areas  of  interest  were 
examined  by  the  model  with  the  following  results. 

Results 

The  results  of  this  thesis  will  be  presented  in  three 
separate  but  highly  related  methods.  The  first  method 
examines  the  data  using  a  one  dimensional  type  of  analysis, 
such  as  the  standard  right  hand  side  ranging  for  linear 
programming.  The  second  method  uses  a  two-dimensional 
analysis  by  comparing  and  plotting  combat  power  as  a 
function  of  two  of  the  independent  variables.  The  third 
analysis  integrates  the  first  two  methods  by  allowing  the 
reader  the  opportunity  to  see  some  of  the  interrelated 
variables  in  terms  of  three-dimensional  pictures. 
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One  Dimensional  Analysis.  APOD  MHE  and  to  a  lesser  extent 
FOL  MHE  were  binding  constraints  on  a  majority  of  the  57 
program  runs,  highlighting  the  importance  of  adequate 
materials  handling  equipment  in  meeting  the  needs  of  the 
theater  commanders . 

Another  important  factor,  which  will  be  expanded  in 
the  two-dimensional  and  three-dimensional  analysis,  was  the 
critical  impact  the  distance  from  the  APOD  to  the  FOL  had 
upon  the  value  of  the  objective  function.  This  point  was 
especially  true  if  the  scenario  did  not  include  any  C-17 
aircraft  in  the  airlift.  Combat  power  decreased  54  percent 
in  one  comparison  when  the  distance  increased  from  zero  to 
10  days  travel  time.  With  40  C-17’s  included,  the  decrease 
was  reduced  to  25  percent.  While  still  a  significant 
figure,  the  findings  highlight  the  negative  impact  of 
distance  and  the  positive  impact  of  C-17  availability  on 
the  delivery  of  combat  power. 

The  number  of  C-130’s  available  was  generally  not  a 
binding  constraint  in  any  of  the  generated  scenarios.  This 
point  was  due  in  part  by  the  limited  ability  of  the  C-130 
to  airlift  a  large  percentage  of  the  cargo  in  this 
particular  scenario.  Over  83  percent  of  the  cargo  by 
weight  consisted  of  either  oversized  or  outsized  cargo, 
neither  of  which  can  be  loaded  on  the  C-130. 

The  C-141  was  utilized  heavily  in  all  the  generated 
scenarios.  Due  to  its  intertheater  and  intratheater 
capabilities,  the  C-141  usage  rates  and  maximum  sorties 

were  almost  always  binding  constraints. 
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The  impact  of  no  fuel  inerting  system  on  the  C-141  was 
examined  from  a  users  standpoint.  The  linear  program  had 
the  advantage  of  knowing  what  the  attrition  rates  would  be 
for  a  specific  scenario  and  the  model  would  adjust  C-141 
usage  accordingly.  For  a  five  percent  scenario  attrition 
rate,  the  model  would  use  10  percent  less  C-141’ s  due  to 
the  expected  loss  rates.  With  a  10  percent  attrition  rate, 
there  would  be  a  44  percent  drop  in  the  C-141  sorties  as 
the  model  attempted  to  minimize  the  loss  of  combat  power 
due  to  aircraft  attrition.  In  reality  attrition  is  not  a 
known  factor,  and  its  impact  upon  the  airlift  force  is  not 
determined  until  the  aircraft  do  or  do  not  return  from 
their  assigned  mission.  Nevertheless,  as  a  point  of 
reference,  if  the  current  aircraft  were  allowed  to  fly  each 
of  the  four  types  of  missions  every  day  of  the  20  day 
operation,  there  would  be  100  percent  attrition  of  the 
C-141 ’s  before  the  end  of  the  operation,  an  unacceptable 
loss  in  terms  of  men,  aircraft,  and  cargo.  A  C-141  flying 
an  airdrop  mission  already  has  the  highest  combined 
attrition  factor,  so  if  the  C-141  is  to  be  used  in  a  role 
other  than  as  an  intertheater  airlifter  from  CONUS  to  APOD, 
a  fuel  inerting  system  should  become  a  high  priority  for 
MAC. 

Two-dimensional  Analysis.  Plotting  combat  power  as  a 
function  of  the  number  of  C-17’s  and  the  distance  between 


the  APOD  and  the  FOL  in  periods  (one  period=5  days), 
revealed  the  importance  of  the  C-17  as  the  distance 


between  the  two  bases  increases.  Setting  the  MHE 
capability  at  1000  pallet  equivalents  per  day  and  aircraft 
attrition  at  5  percent,  figure  five  shows  how  increasing 
the  number  of  C-17's  increases  combat  power  delivered. 

While  zero  days  travel  between  the  APOD  and  FOL  produces 
the  greatest  overall  combat  power  delivered,  the  marginal 
advantage  of  each  additional  C-17  is  highest  when  the 
distance  is  the  greatest,  as  evidenced  by  the  slope  of  each 
line.  With  a  distance  of  zero,  the  marginal  utility  of 
additional  aircraft  is  less,  and  eventually  a  saturation 
point  is  reached  where  other  factors  such  as  MHE  or  ramp 
space  cause  a  decrease  in  combat  power  delivered.  Combat 
power  is  a  function  of  time;  so  if  aircraft  loads  are 
^delayed  due  to  ramp  or  MHE  saturation  the  combat  value  of 
that  load  will  be  reduced.  For  this  particular  scenario, 
with  zero  distance,  the  saturation  point  appeared  to  be 
between  30  and  35  aircraft.  This  saturation  point  would 
occur  with  the  other  two  distances  also,  but  at  a  higher 
level  of  total  aircraft.  The  major  point  of  this  chart  is 
to  emphasize  the  increasing  utility  of  the  C-17  as  the 
intratheater  distances  increase. 

The  plot  of  combat  power  as  a  function  of  MHE  and  the 
number  of  C-17’s  is  depicted  in  figure  six.  As  expected, 
as  the  amount  of  MHE  and  C-17’s  are  increased,  the  amount 
of  combat  power  increased.  The  marginal  utility  of  each 
additional  C-17  increased  at  a  decreasing  rate  (decreasing 
slope)  until  eventually  MHE  capability  is  saturated  and 
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additional  aircraft  provide  no  increased  utility.  The 
tradeoff  between  MHE  and  increased  combat  capability 
appears  linear  throughout  the  range  of  this  graph.  An 
increase  in  MHE  capability  of  300  pallets  per  day  results 
in  an  incremental  increase  in  combat  power  of  12  to  13 
units,  or  approximately  a  4.4  percent  increase  in  terms  of 
the  maximum  amount  of  combat  power  delivered  for  any  one 
scenario . 

A  similar  plot  of  combat  power  as  a  function  of  the 
number  of  C-17’s  and  aircraft  attrition  appears  in  figure 
seven.  As  expected,  the  greater  the  attrition  level,  the 
less  the  incremental  advantage  of  each  additional  C-17. 

The  zero  attrition  level  produced  the  largest  overall 
values  of  combat  power  delivered  as  well  as  the  largest 
marginal  improvement  of  combat  power  per  aircraft. 

As  noted  earlier,  the  C-5  variable  was  not  selected  as 
a  main  effect  in  determining  the  value  of  the  objective 
function.  This  finding  indicates  that  increasing  the 
number  of  C-5’s  from  20  to  60  will  have  no  major  effect 
upon  the  total  value  of  combat  power  delivered.  Using  proc 
stepwise  with  the  "maxr"  option  a  regression  equation  with 
the  C-5  variable  as  a  factor  in  the  model  was  selected  to 
produce  a  regression  equation  where  C-5’s  could  be  compared 
to  distance  (Appendix  D) .  Using  a  response  surface  equation 
with  19  variables,  a  plot  of  combat  power  as  a  function  of 
the  number  of  C-5’s  available  and  the  distance  from  the 
APOD  to  FOL  (figure  8)  verifies  the  earlier  observation 
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that  increasing  the  number  of  C-5’s  beyond  the  20  in  the 
model  produces  negligible  improvement.  While  the  current 
fleet  of  aircraft,  including  the  20  C-5’s,  do  deliver  a 
significant  amount  of  cargo,  the  inclusion  of  up  to  60 
C-5’s  had  a  marginally  positive  impact  on  combat  power  if 
the  distance  from  APOD  to  FOL  is  short.  When  compared  with 
the  C-17  versus  distance  graph  (Fig.  5),  it  can  be  shown 
that  while  the  C-5  delivered  more  cargo,  the  value  of  the 
C-5  decreased  as  the  distance  between  the  bases  increased, 
while  the  value  of  the  C-17  continued  to  increase. 
Three-Dimensional  Analysis.  Several  three-dimensional 
plots  were  produced  using  SAS/Graph,  a  Tektronics  4054,  and 
an  HP7220  plotter.  In  the  first  two  cases  combat  power  was 
examined  while  all  but  two  of  the  main  effects  in  the  model 
were  set  to  their  mid-range  value*  for  the  experimental 
design.  The  results  highlight  several  key  points  discussed 
in  the  previous  analysis.  In  order  to  visualize  the  impact 
of  the  independent  variables,  several  pictures  of  the  same 
plot  were  made  from  different  angles.  Figures  9,  10,  11, 
and  12  depict  the  effects  of  aircraft  attrition  and  the 
number  of  C-17’s  on  the  amount  of  combat  power  delivered. 
The  sequence  of  pictures  rotate  the  plot  approximately  45 
degrees  and  increased  the  tilt  approximately  40  degrees. 

The  maximum  point  in  the  plane  was  the  intersection  where 
C-17’s  were  at  their  maximum  and  attrition  was  at  a  zero 
level.  The  large  decrease  in  combat  power  delivered  can  be 
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visualized  by  observing  the  maximum  point  (237.4),  versus 
the  point  where  attrition  is  10  percent  and  there  are  no 
C-17’s  involved  in  the  operation  (116).  The  decrease 
represented  a  51  percent  reduction  in  combat  power 
delivered  due  to  the  combination  of  high  attrition  and  no 
C-17  aircraft  in  the  scenario. 

The  second  set  of  three-dimensional  plots  was  used  to 
analyze  the  effects  of  distance  and  attrition  on  combat 
power  delivered.  Again  setting  the  other  main  effects  at 
their  mid-range  values,  the  plane  is  depicted  in  figures 
13,  14,  15,  16,  and  17.  The  sequence  rotated  the  plot  45 
degrees  and  left  the  tilt  set  at  approximately  70  degrees. 
The  high  point  of  the  plane  lies  where  distance  equals  zero 
and  attrition  equals  zero.  As  distance  and  attrition 
increased,  the  plane  slopes  downward  and  towards  its  low 
point  where  attrition  and  distance  are  at  the  greatest 
value. 

The  third  set  of  three-dimensional  plots,  figures  18, 

19,  and  20,  depict  the  effect  of  attrition  and  the  APOD  to 
FOL  distance  upon  the  number  of  C-17’s  required  to  reach  a 
certain  combat  power  level.  The  combat  power  level  was  set 
at  a  relatively  low  value  of  145  to  illustrate  the 
incremental  increase  in  requirements  for  the  C-17  as 
attrition  and  distance  were  both  increased  to  ten  percent 
and  10  days  travel  time  respectively.  With  the  other 
aircraft  set  at  their  mid-range  values  of  40  C-5’s  and  60  C- 
130’ s,  the  current  force  was  initially  capable  of  meeting 
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■the  combat  power  requirements.  As  attrition  and  distance 
increased,  the  need  for  the  C-17  became  more  apparent  until 
finally  approximately  35  C-17’s  were  required  to  meet  the 
established  combat  power  value.  These  plots  again 
highlight  the  importance  of  including  the  C-17  in  an 
optimum  airlift  force  structure.  As  the  negative  factors 
of  aircraft  attrition  and  APOD  to  FOL  distance  become  more 
paramount,  so  does  the  requirement  to  include  C-17’s  into 
the  force  structure. 

These  three  sets  of  three-dimensional  plots  not  only 
allow  the  visualization  of  interrelationships  established 
with  the  other  previously  discussed  methods  of  analysis, 
but  they  also  provide  a  means  of  visualizing  the  tradeoffs 
involved  in  determining  a  proper  force  structure  for  a 
given  scenario. 

The  response  surface  equations  used  for  the  analysis 
are  listed  in  Appendix  D.  The  first  equation  is  the  coded 
response  surface  equation.  Any  inputs  to  this  equation 
would  have  to  be  coded  using  the  "+1",  "-1",  and  "0" 
format.  Referencing  the  ranges  applied  to  the  variables 
(see  Table  XIII),  a  specific  input  would  first  have  to  be 
coded.  As  an  example,  30  C-5  aircraft  would  be  coded  as 
( 30-40 )/20  =  -0.5,  with  40  equal  to  the  center  value  and 

20  equal  to  the  difference  between  variables  for  the 
design.  This  coding  method  explains  why  60  C-5’s  was  coded 
as  a  +1  ( 60-40 )/20  =  +1.  The  second  equation  is  the 
decoded  version  of  the  first  equation  so  actual  values, 


within  the  applicable  ranges,  can  be  inserted  into  the 
equation  without  the  requirement  to  code  the  input  first. 
The  third  equation  is  the  coded  response  equation  which 
includes  the  C-5  main  effect  for  comparison  purposes. 
Response.  Surface  Limitations.. 

There  are  several  limitations  to  the  effective  use  of 
these  response  surface  equations,  and  they  include  the 
requirement  to  stay  within  the  bounds  of  the  variables  used 
to  generate  the  response  equation  itself.  Extrapolation 
outside  the  ranges  of  the  factor  levels  can  provide 
erroneous  outputs.  This  first  limitation  highlights  a 
second  restriction.  If  one  of  the  factor  levels  does 
change,  a  new  response  surface  must  be  generated.  A  third 
limitation  is  the  potential  for  requiring  a  large  number  of 
model  runs  to  build  the  response  surface  if  a  large  number 
of  factors  are  to  be  examined  (29:62-63).  A  fourth 
limitation  is  the  hazard  of  bringing  in  a  relatively  minor 
factor  into  the  model  and  using  that  factor’s  minor  effect 
as  a  measurement  of  the  output  of  the  entire  model.  Figure 
eight  demonstrates  this  hazard  by  using  the  minor  main 
effect  of  C-5  aircraft  in  the  model  to  measure  combat  power 
delivered  as  distance  between  APOD  and  FOL  increased.  The 
third  plot  appears  to  show  a  decrease  in  the  objective 
function  value  due  to  the  combined  effects  of  distance  and 
C-5’s.  Yet,  a  linear  program  would  not  show  a  decrease  in 
objective  function  value  as  the  right  hand  side  of  a  less 
than  or  equal  to  constraint  is  increased.  When  viewed  in 


the  context  of  a  maximization  problem  using  linear 
programming,  as  the  amount  of  a  resource  is  increased,  the 
objective  function  will  either  increase  or  remain  the  same 
with  a  corresponding  increase  in  the  slack  value  associated 
with  that  specific  resource.  The  decreasing  objective 
value  was  probably  due  to  a  combination  of  a  less  than 
perfectly  fitting  response  surface  and  the  technique  of 
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Summary 

This  thesis  has  further  refined  a  methodology  which 
combines  response  surface  methodology  and  experimental 
design  concepts  to  describe  the  output  of  a  deterministic 
model.  In  an  effort  to  maximize  the  combat  power  delivered 
to  a  theater  commander,  the  interactive  effects  of  time, 

MHE  capability,  airfield  ramp  space,  distance  between  APOD 
and  FOL,  availability  and  capabilities  of  different 
aircraft,  and  the  impact  of  aircraft  attrition  were 
examined  in  the  model.  Utilizing  an  appropriate 
experimental  design,  a  response  surface  with  a  greater  than 
96  percent  level  of  accuracy  was  generated  for  analysis 
purposes.  The  analysis  demonstrated  the  numerous 
advantages  and  also  the  limitations  of  using  response 
surface  methodology  to  generate  a  multidimensional  analysis 
of  a  given  model.  Even  though  an  optimum  airlift  force  mix 
may  be  scenario  dependent,  this  research  provides 
additional  insight  into  the  complex  interactive  factors 
which  affect  the  ability  of  an  airlift  force  to  meet 
national  objectives. 

CanfitMaigns 

Within  the  ranges  specified  for  this  model,  the 
addition  of  more  than  20  C-5B  aircraft  was  ineffective  at 


increasing  the  combat  power  delivered  to  the  objective 


area.  This  result  does  not  imply  that  the  C-5B  is  not  an 
important  asset,  but  other  variables  such  as  airfield  ramp 
capacity  and  MHE  saturation  prevented  additional  C-5B 
aircraft  from  becoming  an  important  factor.  For  this 
particular  scenario,  with  the  APOD  MHE  limitations  and  ramp 
saturation,  additional  C-5  aircraft  did  not  increase  the 
overall  effectiveness  of  the  airlift  force. 

While  the  location  of  the  battle  may  dictate  which 
airfield(s)  are  designated  as  the  FOL  and  APOD,  the 
distance  between  these  two  bases  was  proven  to  be  a  major 
factor  in  determining  the  quantity  of  combat  power 
delivered  to  the  front,  especially  if  the  C-17  was  not 
included  in  the  airlift  force  structure.  Minimizing  the 
distance  between  the  two  fields  can  significantly  affect 
the  timely  delivery  of  combat  power,  especially  if  units 
must  proceed  overland  from  the  APOD  to  the  FOL. 

APOD  materials  handling  equipment  capability  was  also 
shown  to  be  a  dominant  factor  in  effecting  the  delivery  of 
combat  power.  The  constraint  for  APOD  MHE  was  consistently 
binding  in  the  linear  program,  and  the  response  surface 
equation  shows  APOD  MHE  capability  as  a  primary  factor  in 
determining  the  shape  of  the  response  surface.  Time 
permitting,  the  prepositioning  of  MHE  capability  in  the 
same  theater  as  an  expected  area  of  conflict  could  prove  to 
be  a  major  factor  in  expediting  the  flow  of  men  and 
equipment  to  the  objective  area. 
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The  C-141B  was  shown  to  be  a  valuable  asset  throughout 
the  range  of  the  model.  C-141B  usage  rates  and  sorties 
generated  demonstrated  the  high  utility  of  the  aircraft, 
yet  the  increasing  levels  of  attrition  also  demonstrated 
the  excessive  vulnerability  of  the  aircraft  in  a  threat 
environment.  The  combination  of  using  a  vulnerable 
aircraft  to  fly  on  the  mission  most  susceptible  to  the 
given  threat  (airdrop  mission)  produces  an  extremely  low 
level  of  survivability.  If  the  C-141B  is  to  be  utilized  in 
a  threat  environment,  enhancements  to  decrease  its 
vulnerability,  especially  in  the  area  of  fuel 
vulnerability,  must  be  addressed  and  resolved. 

Aircraft  attrition  was  demonstrated  to  be  a  major 
factor  in  determining  the  quantity  of  combat  power 
delivered  to  the  objective  area.  Combat  power  delivered 
over  the  20  day  operation  averaged  more  than  a  30  percent 
decline  when  the  attrition  factor  increased  to  the  ten 
percent  level.  Even  a  five  percent  level  of  attrition 
produced  decreases  in  combat  power  of  over  20  percent .  The 
response  surface  verified  the  critical  impact  aircraft 
attrition  can  have  upon  the  timely  delivery  of  combat 
power.  While  the  Air  Force  may  not  be  able  to  dictate  the 
susceptibility  of  the  threat  environment,  it  can  increase 
the  aircrafts’  probability  of  survival  by  incorporating 
measures  to  reduce  aircraft  vulnerability  to  the  threat. 


The  C-17  was  shown  to  be  a  significant  factor 
throughout  the  range  of  the  scenario.  Because  of  the 
aircraft’s  flexibility  in  terms  of  delivery  capability, 
inclusion  of  the  C-17  into  the  airlift  force  structure 
provided  an  added  dimension  to  the  capabilities  of  the 
total  force.  While  the  addition  of  40  C-5B  aircraft  did 
not  appreciably  affect  the  quantity  of  combat  power 
delivered  to  the  objective  area  (Figure  8),  a  similar 
addition  of  40  C-17  aircraft  provided  more  than  a  28 
percent  increase  in  the  delivery  of  combat  power  (Figure 
6.)  The  aircraft’s  ability  to  perform  the  airdrop  mission 
combined  with  it’s  capability  to  deliver  outsized  cargo 
directly  to  the  objective  area,  were  major  enhancements  to 
the  current  airlift  force.  As  distance  between  the  APOD 
and  FOL  increased,  the  marginal  utility  of  each  C-17 
increased.  APOD  MHE  saturation  was  also  circumvented  by 
the  direct  delivery  concept  associated  with  the  C-17.  The 
addition  of  this  aircraft  to  the  airlift  inventory  would 
provide  added  flexibility  to  the  operational  planners’ 
capability  to  meet  the  unknown  challenges  of  future 
contingency  operations.  The  impact  of  unknown  factors,  sue 
as  the  APOD  MHE  capability  and  the  APOD  to  FOL  distance, 
could  be  minimized  by  incorporating  the  inherent 
flexibility  of  the  C-17. 


Although  attrition  was  included  in  the  model,  the 
impact  of  aircraft  attrition  upon  the  timely  delivery  of 
combat  power  should  be  examined  in  further  detail.  Rather 
than  using  an  estimation  of  the  combat  power  per  aircraft 
type,  a  more  precise  formulation  should  be  developed  which 
examines  combat  power  by  aircraft  type,  mission  type,  and 
the  period  of  deployment.  Specific  threat  scenarios  should 
also  be  included  for  analysis. 

The  variety  of  units  to  be  deployed  could  be  expanded 
to  include  the  Army’s  High  Technology  Test  Bed  (HTTB)  or  a 
squadron  of  F-15  fighter  aircraft. 

As  stated  previously,  an  appropriate  airlift  force  mix 
is  scenario  dependent.  With  the  current  political 
environment,  consideration  should  be  given  to  developing  a 
scenario  centered  in  South  America.  The  generation  of  a 
new  scenario  could  be  expedited  by  the  development  of  a 
matrix  generator  to  build  the  linear  model. 

The  impact  of  the  C-130  appears  to  have  been  decreased 
by  the  high  percentage  of  oversized  and  outsized  cargo  in 
this  particular  scenario,  but  further  analysis  is  required 
to  verify  that  conclusion. 

While  the  current  analysis  highlighted  the  critical 
impact  of  the  distance  between  the  APOD  and  FOL,  the  impact 
of  the  distance  from  the  Aerial  Port  of  Embarkation  (APOE) 
and  the  APOD  could  also  be  an  important  factor  and  may  be  a 
source  for  further  analysis. 
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Appendix  C 


SAS  REGRESSION  PROGRAM 


option*  linesize  ■  76; 
filenaae  newdata  'design.dat') 
data  first; 

infila  newdata; 
labtl  x 1 - ' C-5  AIRCRAFT'; 
labal  x2«'C-17  AIRCRAFT’) 
label  x3* ' APOD  MHE't 
label  x4«'ilRCRAFT  ATTRITION') 
label  x3* ' DIST  TO  FOL'i 
label  /6-'AP0D  RAMP  SIZE') 
label  x7* 'C-130  AIRCRAFT') 
label  y* ' COMBAT  FIRE  POWER'; 
input  xi  x2  x3  x4  x5  x6  x7  y  j 


x 1 *x 1 ) 

■  xl*x2) 

»  xl*x3; 

■  x  1  *x  4  j 

■  x I*x5j 

■  xt#x6j 

*  xl#x7j 

■  x2*x2j 

■  x2*x3) 

■  x2#x4j 

■  x2*x5) 

■  x2*x6) 

*  x2*x7i 

■  x3*x3 j 

■  x  3  *x  4  j 

*  x3*x5) 

■  x  3  •  x  6 ; 

■  x3*x7j 

■  x4*x4) 

■  x4*x5j 

■  x4*x6) 

■  x4*x7) 

■  x5*x3) 

■  x3*x6) 

■  x5*x7; 

■  x6*x6) 

■  x6»x7j 

■  x7*x7; 
proc  stepwise; 

aodel  y  «  xl  x2  x3  x4  xS  x6  x7  xii  x  1 2 

x 13  x 14  x 13  x 16  x 17  x22  x23  x24 

x23  x 26  x 27  x33  x34  x35  x36  x37 

x44  x43  x46  x47  x53  x56  x37  x66  x67 

x77/taxr ) 


xll 
x  12 
x  13 
x  14 
xl5 
x  16 
xl7 
x22 
x23 
x24 
x23 
x26 
x27 
x  33 
x34 
x33 
x36 
x37 
x44 
x45 
x46 
x47 
x53 
x36 
x57 
x66 
x67 
x77 


/»  Regress  on  ALL  Variables  * 

proc  reg  data  ■  first) 

Model  y  ■  xl  x2  x3  x4  xS  x6  x7  xll  x  1 2 

x 13  x 14  x IS  x 16  x 17  x22  x23  x24 

x23  x26  x27  x33  x34  x35  x36  x37 

x44  x45  x46  x47  x33  x36  x37  x66  x67  x 77 ; 

output  out«resids  p«pred  r«resid) 

data  residul) 
set  resids) 
errorl  ■  resid/pred) 
err  1  ■  abs(errorl)! 
proc  plot  data  *  residul) 

plot  resid*pred  »  "RM/vref «0j 
proc  plot  data  ■  residuii 

plot  errorl*pred«  r,EM/vref«  Oj 
proc  print; 

var  errorl  y  pred) 
proc  eeans  Mean  ain  nax] 
var  y  pred  errl) 


/* 


proc  reg 
aodel  y 


data 


output  out  *  resoids2 
data  residual 


Regress  on  Selected  Model 
*  firsti 

x2  x3  x4  x3  x 15  x22  x24  x25 
x34  x33  x44  x53/corrbj 


ta  residual) 
set  resoids2) 
error  ■  res/prei 
err  ■  abs(error)) 
proc  plot  data  *  residual! 

plot  error*pre/vret«Oj 
proc  print) 


pre  r  ■  resj 


var  error  x2  x3  x4  xS  x 1 3  x22  x24 

x23  x34  x35  x44  x35  y  pre) 
proc  Means  aean  ain  aax  ; 

var  x2  x3  x4  x5  xl5  x22  x 24  x23 
x34  x33  x44  x33  y  pre  err) 


RESPONSE  SURFACE  EQUATIONS 


Cod»d  BlQflMi  Surface  Equation  (12  Variable)  i 


y  -  189.986739  +  22. 250625 X_  +  12. 78904 IX,  -  3B.021625X. 

2  3  4 

-  43.942458XS  -  5.93073X^3  -  12.717293X22  -  6.713X2X4 
♦  12.361X2Xs  -  3.71075X3X4  -  8.7935X3XS  +  5.026456X44 
«•  23. 329331 X55 


Ctanyr  tad  BMBflDM  Surf  ace  Eauat  ion  (12  Var  tables)  i 


y  -  174.97222  +  2. 101961X2  +  .084311X3  -  379.80743X4 

-  62.14B955XS  -  .29753BX1XS  -  .031793X22  -  6.713X2X4 

+  .61803X2X3  -  .247383X3X4  -  . 02931 167X3XS 

+  2010.382608X  2  +  23.329332X  2 
4  3 


Coded  BlBflOll  Surface  Equation 


<To  include  C-5  win 


iff  met.  19  variables) » 


y  m  194.473134  +  .779541XJ  +  22.230625X2  +  12.789041X3 

-  38.021625X4  -  43.942438XS  +  1.234373X6  -  3.130048XJ1 

-  2.447XjX3  -  5.95075X^3  -  14.119298X22  +  2.335625X2X 

-  6.713X2X4  +  12.361X2X5  -  3.299173X33  -  3.71075X3X4 

-  8.7933X3Xs  +  1.369875X3X6  3.624451X  4 
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